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TORIC NETWORKS, GEOMETRIC R-MATRICES AND 
GENERALIZED DISCRETE TODA LATTICES 

REIINOUE, THOMAS LAM, AND PAVLO PYLYAVSKYY 


Abstract. We use the combinatorics of toric networks and the double affine geometric 
i?-matrix to define a three-parameter family of generalizations of the discrete Toda 
lattice. We construct the integrals of motion and a spectral map for this system. 
The family of commuting time evolutions arising from the action of the i?-matrix is 
explicitly linearized on the Jacobian of the spectral curve. The solution to the initial 
value problem is constructed using Riemann theta functions. 


1. Introduction 

We study a generalization of the discrete Toda lattice parametrized by a triple of 
integers (n, m, k), which corresponds to a network on a torus with n horizontal wires, m 
vertical wires and k shifts at the horizontal boundary. An example of the kind of toric 
network that we consider is given in Figured! 



Figure 1. The toric network for n = 3, m = 4 and /c = 0. 

The phase space M. of our system is the space of parameters Qij G C assigned to 
each of the crossings of the wires (a 3 x 4 = 12 dimensional space for the example in 
Figured])- We construct two families of commuting discrete time evolutions acting on 
the parameters qij, together generating an action of Z™ x where N := gcd{n,k). 
These time evolutions act as birational transformations of the phase space. The purpose 
of this paper is to study the algebro-geometrical structure of these maps and to solve 
the corresponding initial value problem. 

The rational transformations of the phase space come from the affine geometric R- 
matrix, acting on the parameters of two adjacent parallel wires, either horizontal or 
vertical. The affine geometric R-matrix arises in the theory of affine geometric crystals 
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[BKllKNT|IKN2j . being a birational lift of the combinatorial R-matrix of certain tensor 
products of Kirillov-Reshetikhin crystals for 17q(sR) [KKMMNNj . Geometric /^-matrices 
also arises independently in the study of Painleve equations |KNYj and total positivity 
pVPT^ : see also jKiFllEti] . 

The geometric i?-matrix satishes the Yang-Baxter relation and we show in Theorem 
12.41 that the action of the geometric i?-matrix generates commuting birational actions 
of two affine Weyl groups W and W\ one swapping horizontal wires, and one swapping 
vertical wires. The commutativity was proved in |KNYj for the case A; = 0; we extend it 
here to arbitrary k. 

The X time-evolutions of our dynamical system come from the subgroup of 
translation elements in the corresponding affine Weyl group. In the case (n, m, k) = 
{n,2,n — 1), a rational map given by one of the Z™^^ actions corresponds to the dis¬ 
cretization of the well-known n-periodic Toda lattice equation, studied in |HTIj . 

The dynamics that we study come from positive birational maps, and they tropicalize 
to the box-ball system [TSj and the combinatorics of jeu-de-taquin (see ^2.5.2p . 

To study the dynamics of our generalized discrete Toda lattice, we construct a spectral 
map 

4 >: A4 —;■ {plane algebraic curve} x Pic^(G/) x Sf x Rq x Ra 
(qij) I— {Cf, V, (ci, ...,cm),0,A) 

where G/ is a spectral curve, V G Pic^(G/) is degree g divisor, and the remaining data 
is explained in ^ We show that when appropriately restricted, the spectral map is 
an injection. Such spectral data is frequently encountered in the theory of integrable 
systems, and our approach follows that of van Moerbeke and Mumford |vMM] . Van 
Moerbeke and Mumford used similar spectral data to study periodic difference operators. 

The heart of our work is the calculation of the double affine geometric i?-matrix action 
in terms of the spectral data: we show that the translation subgroup Z"* x Z^ acts as 
constant motions on the Jacobian of Cf, while the symmetric subgroups of W and of W 
act by permuting the additional data Rq and Ra (which are certain special points on 
Cf). 

We explicitly invert (for N = 1) the spectral map 0 using Riemann theta functions, 
and give a solution to the initial value problem. This extends work of Iwao |I08llI10j . 
who studied the initial value problem of the Z”^ action in the (n, m, n — 1) and {n, m, 0) 
cases. We relate our theta function solutions to the octahedron recurrence |Spe| via 
Fay’s trisecant identity 1^ - We also give an interpretation of our dynamics in terms 
of dimer model transformations EKl. 

Outline. The outline of this paper is as follows: in l|2l we introduce a family of com¬ 
muting actions on the toric network, generalizating the discrete Toda lattice. We briefly 
summarize the main results of this paper in 1 12.41 and give some examples and applica¬ 
tions in 1 12.51 In 1|31 we introduce a space C of Lax matrices for our system, and construct 
the spectral curves Cf, where / G C[x,y]. We also study certain special points on our 
spectral curves, and analyze some of the singular points. In 1|51 we study the spectral 
map which sends a point in the phase space to the spectral curve Cf, a divisor V on 
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Cf, and some additional data. We modify the strategy in |vMMj to fit our situation. In 
§4] and §6] we present the proofs of results in ^ and ^ respectively. In particular, the 
coefficients of the spectral curves (which are integrals of motion) are described explicitly 
in terms of the combinatorics of these networks. In we study the vertical actions, 
the horizontal actions, and the snake path actions. We show that all the actions 
preserve the spectral curve, and the commuting m + N actions induce constant motions 
on the Picard group of Cf, through the map 0. In ^ we solve, for the case of = 1, 
the initial value problem for the commuting actions by constructing the inverse of 0 in 
terms of the Riemann theta function. Our method extends the strategy in |Iin] . For 
TV > 1, a solution is given that relies on a technical condition. In §21 we show that the 
theta function solution for the system satishes the octahedron recurrence, by specializing 
Fay’s trisecant identity for the Riemann theta function. In I JTOl we study the symmetry 
between the network and its transposition obtained by swapping the roles of the vertical 
and the horizontal wires. In I fTTl we realize the R-matrix transformation on a toric net¬ 
work in terms of transformations of the honeycomb dimer model on a torus EE]. An 
explicit interpretation of the R-matrix as a cluster transformation will be the subject of 
future work ra. so we do not elaborate on the relation to cluster structures here. 

Future directions. There are several systematic ways to construct integrable ratio¬ 
nal maps using combinatorial objects on surfaces, such as directed networks, electrical 
networks, and bipartite graphs |LP131 ITjK 1IGSTV] . In particular, the R-matrix of the 
present work has an “electrical” analogue |LP12b] . It would be interesting to study these 
discrete-time dynamical systems from the view point of algebraic geometry (Cf. [FMj j. 
In all these cases, the rational transformations are additionally positive, and can be 
tropicalized to piecewise-linear maps that generate a discrete-time and discrete-space 
dynamical system. It is natural to consider the tropical counterpart of our work using 
tropical geometry as in nn. where tropical curves and tropical theta functions are used 
to study the piecewise-linear map arising from the {n, m, k) = {n, 2,n — 1) case. 

Acknowledgments. We thank Rick Kenyon for helpful discussions. We also thank the 
anonymous referee for kind comments which improved the manuscript. 


2. Dynamical system from a toric network 

We study the action of the geometric R-matrix on an array of variables. The geometric 
R-matrix generates an action of a product of commuting (extended) affine symmetric 
groups, acting as birational transformations. 

2.1. The afRne geometric R-matrix. For a vector a = (ai,...,a„), let a^ ; = 
(oj+i, ai+ 2 ,..., a„, oi,..., Oj). Let a = (oi,..., a„) and b = (6i,..., 6„) be two vectors. 
Dehne the energy E{sl, b) to be 

n—1 / i n 

E{a,h )=XI n 

i=0 \j=l j=i+2 
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Define the affine geometric i?-matrix to be the transformation R : (a, b) i—)■ (b', a'), given 
by 


, _ h;(a«,b«) 

* “ ^*E(a(*-i),bh-b) 


and 


E(aW,bh)) 


We will usually think of i? as a birational transformation of x C"'. 
It is easy to see 


n =n n =n 

i=l 2=1 2=1 2 = 1 


2.2. Description of dynamics. The i?-matrix can be used to act on a rectangular 
array of variables, by acting on consecutive pairs of rows or of columns. 

Let a, b, c be three positive integers, and let 0 < c) < c be a nonnegative integer 
satisfying gcd(c, b) = 1. Note that we allow b = 0 if and only if c = 1. Let b“^ be the 
unique number in the range 0 < b“^ < c such that bb”^ = 1 mod c. For a positive 
integer a, we define [a] := {1, 2,..., a}. 

We consider an array of ax b x c variables {(?a,fe,c}ae[a],fee[6],ce[c]- We consider two distinct 
ways to arrange them in a two-dimensional rectangular array: 

Qij = qa,b,c ioT i = a, j = b + b(c - 1) (2.2) 

and 


Qi,j = qa,b,c for i = b -|- 1 — b, j = ab — a -|- 1 mod ac. (2.3) 

By convention, gjj denotes the entry in the i-th column and j-th row. We call these 
arrays Q and Q. The first array Q has dimensions be x a and the second array Q has 
dimensions ac x b. Note that changing a i-a b, b i-A b~^, qa,b,c t Q'i-6,i-a,a-ic swaps Q 
and Q. 


Example 2.1. Let (a, b, c, b) = (2, 2, 2,1). The two arrays are as follows. 


Q: = 


( ?2,1,1 \ 
qi,2,l Q2,2,1 
qi,l,2 Q2,1,2 
\ qi,2,2 ^2,2,2 / 


Q 


^ Q2,2,1 Q'2,1,1 ^ 

Qi, 2,1 qi,i,i 
Q2,2,2 Q2,1,2 
\ <11,2,2 qi,l,2 / 


Example 2.2. Let (a, b, c, b) = (3, 2, 3, 2). The two arrays are as follows. 



'? 2 , 1,1 

^ 3 , 1,1 

\ 

qi , 2,i 

<? 2 , 2,1 

^ 3 , 2,1 


Qi , 1,2 

<? 2 , 1,2 

93 , 1,2 


qi , 2,2 

'? 2 , 2,2 

93 , 2,2 



<? 2 , 1,3 

93 , 1,3 


V ? 1 , 2,3 

? 2 , 2,3 

93 , 2,3 

/ 


/ ^ 3 , 2,1 

93 , 1,1 

\ 

92 , 2,1 

92 , 1,1 


9 i , 2 ,i 

9l , l,l 


93 , 2,3 

93 , 1,3 


92 , 2,3 

92 , 1,3 


9 i , 2,3 

9 i , i ,3 


93 , 2,2 

93 , 1,2 


92 , 2,2 

92 , 1,2 


V '? 1 , 2,2 

9 i , i ,2 

/ 


(2.4) 


(2.5) 


The arrays Q and Q correspond to a network on a torus in the following way. Let G 
a network embedded into a torus, as illustrated in Figure [H Each rectangle of red lines 
denotes the fundamental domain of the torus. There are a vertical “wires” (directed 
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upwards), forming a simple curves in the torus with the same homology class, and be 
horizontal wires (directed to the right), which form b simple curves on the torus with 
another homology class: the top bb right ends of horizontal wires cross the top edge 
and come out on the other side. We set qij at the crossing of the f-th vertical and the 
j-th horizontal wires, then Q corresponds to the conhguration of the qij on G. When 
we see G from the inside of the torus, the roles of the vertical and horizontal wires are 
interchanged; there are b vertical wires forming b simple curves, and ac horizontal wires 
forming a simple curves. The top ends of horizontal wires cross the top edge and 
come out on the other side. The array Q corresponds to this “from-inside” conhguration, 
where is placed at the crossing of the f-th vertical and the j-th horizontal wires. 

See Figure [3] for the network in the case of Example 12.21 where the lines with the same 
color in the two pictures are identical simple curves in G. 
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Figure 2. Toric network 


Definition 2.3. Let A = {aij}i^[s],j£[r] be an r x s array. Let a* = {aij)j^[r], and write 
slJ for the transpose of a*. For 1 < ^ < s — 1, dehne the array S£{A) by 

se{A) = (a[,..., aj_i, , (a^)^, aJ+ 2 , • • •, aj) (2.6) 

ii A = (af,..., aj, aJ_^_-^^, aj_^ 2 ) • • •) af), and R{ai, a^+i) = (a^^;^, a^). Also dehne an 

array 7r(A) by the formula 

~ O'iq-i- (2-7) 

Here we consider the hrst index modulo r and the second index modulo s. 

Let M. ~ be the phase space of our dynamical system and regard (ga,6,c)ae[o],fee[6],ce[c] 
as coordinates on M.. Applying the above dehnition to the array Q, we obtain operators 
si,..., Sa_i and tt acting on M.. Similarly, applying this dehnition to the array Q, we 
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Figure 3. Network for the case of (3,2,3,2) 


obtain operators si,..., st,_i and tt acting on A4. We emphasize that fl2.7p for Q gives 
the operation tt, while fl2.7l] for Q gives the operation vr. 

The following resnlt generalizes a result of Kajiwara, Noumi and Yamada [KNY] . 

Theorem 2.4. The operators S£ and vr generate an action of an extended affine symmet¬ 
ric group W = (Z/acZ) ix ©„ on M.. Here 'LlacL is the cyclic group of order ac and is 
the affine symmetric group (the Coxeter group) of type Aa-i. Similarly, the operators §£ 
and 71 form an action of an extended affine symmetric group W = (Z/bcZ) k ©t, on M., 
where (Z/bcZ) is the cyclic group of order be. Furthermore, these two actions commute. 
Precisely, on M the following relations hold. 

SjSi S)/Sj (|j '^1 ^ 1)5 S£ 1, 

TTSi+i = SiTl, 7r“' = 1, 

S£S£-\-iS£ , SjSi S^Sj (|j '^1 ^ 1)) S£ 1, 

TTSf+l = S^TT, = 1, 


SjS£ = SiSj, TlSi = Si'll, 71 Si = SiTl, 7171 = 7171. 

In the above formulae, the index of Si is taken modulo a, and Sq is defined by the eguation 
7isi = Sovr. Similarly, the index of si is taken modulo b, and sq is defined by the eguation 

TTS\ = SqTT, 


The proof of Theorem 12.41 is delayed to 1 14.51 


Example 2.5. In Example 12.11 we have 

+ <?2,2,1<?1,2,25'1,1,1 


'Sl(Q)l,l,l — <?2,1,1 
'Sl(Q)l,l,l = <?1,2,1 


Q'2,2,iQ'2,1,25'1,1,1 + <?2,2,1<?2,1,2<?2,2,2 


Q'1,2,1<?1,1,2<?1,2,2 + <?2,1,1<?1,1,25'1,2,2 + Q'2,l,lQ'2,2,lQ'l,2,2 + <?2,1,1<?2,2,1<?2,1,2 ’ 

5'1,1,1<?2,1,2<?1,1,2 + <?1,1,1<?2,1,25'2,2,1 + 5'l,l,l5'l,2,2Q'2,2,l + <?2,2,2<?1,2,2<?2,2,1 

?2,1,2<?1,1,2<?2,1,1 + <?2,1,2<?1,1,2Q'1,2,1 + Q'2,l,2Q'2,2,lQ'l,2,l + <?1,2,2<?2,2,1<?1,2,1 ’ 
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~ 51,2,1’ ^(Q)l,l,l — 52,1,2- 

We leave it for the reader to verify that si and tt commute with si and tt, as an easy 
computational exercise. 

We can present the extended affine symmetric groups W and W as follows. The hnite 
symmetric group ©o acts naturally on the lattice Z“, hxing the subgroup generated by 
the vector (1,1,..., 1). Thus we have an action of ©„ on the quotient Z“/Z(c, c,..., c). 
Then we have that W = x (Z“/Z(c, c,..., c)). The commutative normal subgroup 
Z“/Z(c, c,..., c) is generated by e„ for 1 < n < a: 

('U (^u ■ ■ ■ -So—1) (-^n—1 ■ ■ ■ - 50 — 2 ) ■ ■ ■ (-^l ■ ■ ■ ^a—u)^ ■ (^•®) 

The element is identihed with the vector e„ = (1,..., 1, 0,..., 0) G Z“ with u 1- 
s. Note that Ca = 7r“ satishes e„ = 1 agreeing with the fact that (c, c,..., c) = 0 in 
Z“/Z(c, c,..., c). Similarly, we have W = ©[,k(Z^/Z(c, c, ..., c)), where the commutative 
normal subgroup Z^/Z(c, c,..., c) is generated by e„ for 1 < u < b: 

Cn (-5^ ■ ■ ■ -Sb—l)('Su—1 ■ ■ ■ -56—2) ‘ ' ' (-Si ' ' ' -^b—• (^•^) 

Now, the operators and all commute, and thus give an action of Z“/Z(c, c,..., c) x 
ZVZ(c,c,...,c) on M.\ We will think of this as a discrete dynamical system with a + b 
different time evolutions. We shall hnd a complete set of integrals of motion for this 
system, and study the initial value problem. 

2.3. Change of indexing. Instead of the quadruple (a, b, c, b), we shall also index our 
systems with triples (n, m, fc), given by 

n = be, m = a, k = bb. 

We shall also set 

N ■= gcd(n, k) = b, n' := n/N = c, k' := k/N = b, M := gcd(n, m + k). 
The quadruple (a, b, c,b) can be recovered via 

a = m, b = gcd(n, fc), c = n/gcd(n,/c), b =/c/gcd(n,/c). 

The involution Q 1 —> Q is associated with the following changes of indices: 

(a, b,c,b) I—^ (b,a, c,b~^), 

(n, m, k) I —> {mn', N, k'm), 
where k' := {k')~^ is taken modulo n'. 

Through 02.21) or 02.31) we identify q = {qa,b,c) ^ -M with Q = (qij) G Matn,m(C) 
or Q = {qij) G Matin', a(C). Correspondingly, the network G has m vertical wires 
(directed upwards) which form m simple curves on the torus with the same homology 
class, and n horizontal wires (directed to the right) which form N simple curves on the 
torus with another homology class. 
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2.4. Main results. Let A4 ~ = C™'"' be the phase space where the ring 0(A4) of 

regular functions on M is generated by qij (i G [m], j G [n]). In ^ we shall dehne a 
map 'ip : Ai ^ £.[x,y\. Every coefficient of f{x,y) = 'ipi.q) is a regular function on M.. 

We then have (see Corollary 17.6p the following result. 

Theorem. The actions of the affine symmetric groups W and W on Ai preserve each 
fiber 'ip~^{f) for f G ^^{Ai). In particular, the coefficients of f{x,y) are integrals of 
motion of the commuting Z™ and Z^ actions. 

In ^3.21 we give a combinatorial description of every coefficient of /, as generating 
functions of path families on a network on the torus. 

Now fix a generic / G 'ip{Ai) such that the affine plane curve {(x, y) \ f{x, y) = 0} C 
is smooth (except for (0,0)), and let Cf be the smooth completion of the affine curve. 

We shall call Cf the spectral curve. In ^ we define distinguished special points, P, 

Au {u G [m]) and Ou {u G [N]) on Cf. We apply the results and techniques of van 
Moerbeke and Mumford |vMM] to establish the following; see Theorem 15.41 

Theorem. Fix a generic f G 'ip{Ai). There is an injection 

cf : fj-\f) -A Pic^(C/) xSfXRoX Ra, 

where 

(1) g is the genus ofCf, and PiC{Cf) is the Picard group ofCf, of degree g; 

(2) Sf ~ (C*)^-i C (C*)^; 

(3) Ro is a finite set of cardinality N\, identified with the N\ orderings of {Oi, O 2 ,... ,On}; 
and 

(4) Ra is a finite set of cardinality m\, identified with the ml orderings of{Ai, A 2 ,..., Am}. 
In fact, our Theorem 15.41 identifies the image of (p. 

Define the divisors Au ■= uP—Yf^=i Ai and Ou '■= uP—'Y}!j=N-u+i ^j- The commuting 
time evolutions e„ fl2.8p and fl2.9p can be described as follows (Theorem 17.21 and 
Theorem 17.3p . 

Theorem. Suppose (j){q) = ([D], (ci,..., cm), O, A). Then we have 

0(en(g)) = {[V- Au], {cu+i, ...,cm,ci,...,Cu), 0,A) for u = 1,..., m, 

4>{eu{q)) = {[D + Ou], {cm-u+ 1 , ..., cm, ci,..., cm-«), 0,A) for u = 1,..., N. 

Furthermore, the finite symmetric subgroups ©at C W and &rn C W act naturally on 
Ro and Ra respectively and do not affect the rest of the spectral data. 

In other words, the time evolutions e„ and e„ are linearized on PiC[Cf). Our approach 
to Theorem 17.31 is similar to that of Iwao |I08j . 

When = 1, we give in Theorem 18.51 an explicit formula for the inverse to the map 
(j), which also explicitly solves the initial value problem for our dynamics. For t G Z™, 
let (gj j) denote the point in the phase space after time evolution in the direction t. 

Theorem. When N = 1, we have a formula 

Ij.i 
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where 6\{P) is a particular value of the Riemann theta function, and C is a constant 
depending only on (ci,..., cm), O, and A. 

2.5. Examples. 


2.5.1. Discrete Toda lattice. Let (a, b, c, O) = (2,1, n, n —1) (i.e. {n,m,k) = {n,2,n — l)). 
We set Qij = qa,i,c for i = a, j = c, and regard q := {qij)i^i 2 ],je[n] as a coordinate of 
M ~ C^"'. The resulting Z^-action on M is generated by ei and 62 , where 62 simply acts 
on M as e 2 {qij) = qij+i- As for ei, when we dehne q^ := e\{q) for g G AT, the action of 
Cl is rewritten as a system of difference equations: 


qpAi + AA = Aj+i + ^Li+i- 


( 2 . 10 ) 


This is the discretization of n-periodic Toda lattice equation studied in |HT1] . Actually, 
we recover the original Toda lattice equation -^Xj = Py setting q\j = 

1 + S-^Xj and q^ j = and taking the limit 5^0. Here we set ql j = q^:j{6t) 

for * = 1, 2 and Xj = Xj(t). 

A simple generalization of the discrete Toda lattice is the case of (a, b, c, b) = (m, 1, n, n— 
1). Its initial value problem was studied by Iwao 110811110] by applying [vMM] . He also 
studied the similar problem in the case of (a, b, c, b) = (m, 1, n, 0) in |109j . 


2.5.2. Tropicalization and tableaux. Let Ti and T 2 be two semistandard tableaux of rect¬ 
angular shapes. Dehne Ti®T 2 to be the concatenation of the two into a skew tableaux by 
placing T 2 North-East of Ti. Let jdt(Ti( 8 )T 2 ) denote the straight shape tableau obtained 
by performing jeu-de-taquin on Ti 0 T 2 . The following result is well-known. 


Lemma 2.6. There is a unique pair of rectangular semistandard tableaux T[ and T 2 such 
that TI has the same shape as R (for i = 1,2), and jdt(Ti 0 T 2 ) = jdt(T 2 0 T[). 


Example 2.7. Suppose 


Ti = 


and 


To = 


12 3 


Then one has 


T' — 
-‘-2 — 


1 3 3 


and 


n 


1 ^ 

2 


2 

3 


The transformation R(Ti,T 2 ) = (T 2 ,T() is called the combinatorial R-matrix. It ap¬ 
pears in the theory of crystal graphs as the isomorphism map between tensor products 
of Kirillov-Reshetikhin crystals. The following property is well-known. 


Theorem 2.8. The combinatorial R-matrix is an involution. Furthermore, it satisfies 
the braid relation: 

{R® Id){Id® R){R® Id) = {Id® R){R® Id){Id® R). 

Now, let {ga,6,c}ae[a],6e[6],ce[i:] be an array of nonnegative integers. Dehne Q = {qij) and 
Q = ilij) before: 


qi,j = qa,b,c ioT i = a, j = b + b(c - 1 ) 
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and 

Qij = (la,b,c for i = b + 1 — 6, i = — a + 1. 

We create single-row tableaux from Q and Q as follows. For each i G [a], let T, be the 
single-row tableau with j-s. For each i G [b], let T* be the single-row tableau with 
Qij j-s. We view T = Ti 0 ... 0 Tq and T = Ti (g)... (g) Tf, as tropical analogues of Q and 

q’ 

Let us dehne an action of W on qa,b,c as follows. For T = Ti (g)... (8)Tn and 1 < £ < a — 1 
let 


se{T) = Ti (g)... (g) 0 T; 0 ... 0 t„. 

That is, we apply the combinatorial i?-matrix to the .^-th and {i + l)-st factors. In 
addition, let tt act on T by applying Schutzenberger’s promotion operator |Sch] to each 
factor Tj. Similarly, we dehne S£ and vr acting on T. 


Theorem 2.9. The operators se and tt form an action of the extended affine symmetric 
group W on qa,b,c- Similarly, operators si and tt form an action of the extended affine 
symmetric group W on qa,b,c- These two actions commute. These two operations are the 


tropicalizations of the rational actions of Theorem \2.4 


Here tropicalization refers to the formal operation of substitution 
C I — y Z, X I —} -j-, -G I —} —, -j- I — y min . 


Under this substitution, the geometric i?-matrix becomes a piecewise-linear involution of 
Z" X Z”, which is the combinatorial i?-matrix. In other words, the W dynamics we are 
considering in this paper is a birational lift of the dynamics of repeated application of 
the combinatorial i?-matrix on a sequence of a single-row tableaux, arranged in a circle. 


Example 2.10. Let (a, b,c,b) = (2,2,2,1) as in Example 12.11 Take T = 22222344 0 
1234. Then T = 122222344 0 134. We have 

si(T) = 2234 0 12222344, si(T) = 111122334 0 134; 
fr(T) = 11111233 0 1234, fr(f) = 123 0 122222344; 
si(T) = 11112334 0 1233, si(T) = 124 0 122223344; 

7r(T) = 1234 0 22222344, 7r(T) = 111112334 0 234. 

Remark 2.11. The corresponding commuting crystal actions in the case c = 1 were 
considered by Lascoux in |Lasj and by Berenstein-Kazhdan in |BKj . 


2.5.3. Box-ball systems. The box-ball system is an integrable cellular automaton intro¬ 
duced by Takahashi and Satsuma na, It is described by an algorithm to move hnitely 
many balls in an inhnite number of boxes aligned on a line, where a consecutive array of 
occupied boxes is regarded as a soliton. This system is related to both of the previous 
two examples; the global movements of solitons are equivalent to the tropicalization of 
the discrete Toda lattice 02.101) . The symmetry of the system is explained by the crystal 
base theory, where the dynamics of bolls is induced by the action of the combinatorial 
R-matrix. 

See |IKTj for a comprehensive review of the combinatorial and tropical aspects of the 
box-ball system. 
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3. Lax matrix and spectral curve 


3.1. Lax matrix. Fix integers n,m,k such that n > 2, m > 1 and 1 < k < n. From 
now on we shall mainly use q := (q’jjOieHjeW as a coordinate of the phase space At. We 
identify g G At with an m-tuple of n by n matrices Q := with a spectral 

parameter x, where 


Qi{.x) 


/ gii 0 0 X \ 

1 gi2 0 0 

0 ■•. ■•. 0 

\ 0 0 1 Qin J 


(3.1) 


Let C be the set of n by cxo scalar matrices A := {aij)i<i<n,jez satisfying the following 
conditions: 


{ 1 j — i = —fn — k 

ttij G C —m — k + l<j — i< —k (3-2) 

0 otherwise. 


In particular, A has hnitely many nonzero entries. For A G £, we dehne an n by n 
matrix L{A;x) = {l{x)ij)i<ij<n by 


eez 


(3.3) 


We may identify C with C™", and identify A E C with L(A; x). 
Now dehne a map a ; At —)■ £ by 


a: Q = {Qi, Q 2 , • • •, Qm) I—^ L{x) := Qi{x)Q 2 {x) ■ ■ ■ Qm{x)P{xf, 

where 



0 

0 

X \ 

1 

0 

0 

0 

0 



0 

V 0 

0 

1 

0 / 


(3.4) 


We call L{x) the Lax matrix. Our approach to the study of L{x) is close to that of van 
Moerbeke and Mumford |vMMj . 

We give a combinatorial description of the Lax matrix by using highway paths on the 
network G. We introduce the “x-line” as illustrated in Figure 01 where the top k right 
ends of horizontal wires cross the x-line (the top edge) and come out on the other side. 
There are n sources labeled 1, 2,..., n on the left and n sinks labeled 1, 2,..., n on the 
right. Each sink j is connected to the source j + k mod n by a wire, as illustrated in the 
hgure. There are mn intersection points between the vertical wires and the horizontal 
wires, which we call the {i, j)-crossroads, where i = 1,2,... ,m indexes the vertical wire, 
and j = l,2,...,n indexes the horizontal wire. Let us denote by G' the cylindrical 
network obtained by gluing only the upper and lower edges of Figure 01 In the torus 
network G, source i and sink i are the same point. In the cylindrical network G', source 
i and sink i are distinct points. 

Let us introduce the notion of highway paths following |LP13j . A highway path p is 
a directed path in the network G (or G') with the following property: at any of the 
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X-line 



Figure 4. The toric network G, with a coil and a snake path shown. 

crossroads, if the path is traveling upwards, it must turn right. We shall only consider 
highway paths that start at one of the sources 1,2,... ,n and end at one of the sinks 
1,2,... ,n. The weight wt(p) of a highway path p is defined as follows. Every time p 
passes the x-line it picks up the weight x. Every time p goes through the (z, j)-crossroad, 
it picks up the weight qij or 1, according to Figure |5l The weight wt(p) is the product 
of all these weights. The condition that a highway path must turn right when travelling 
up into a crossroad is indicated in the Figure [S] we can think of such a turn as giving 
weight 0. Finally, a highway path p may use no edges. In this case, we consider the path 
p to start at some source i, and end at the sink i. We declare such paths to be abrupt, 
and have weight wt(p) = —y. 


weights: 1 q^j 1 0 

Figure 5. Highway paths 

The following lemma gives a highway-path description of the Lax matrix, which is a 
variant of the results of |LP13] . It follows directly from the definitions. 

Lemma 3.1. Let L{x) = a{q) for g G Al. For 1 < i, j < n, we have 

{i,j)-th entry of L{x) - y = ^wt(p), 

p 

where the summation is over highway paths in G' from source i to sink j. 

3.2. Spectral curve and Newton polygon. Dehne a map ^ : A4 — ?■ C[x,y] as the 
composition of a : A4 £ and the map ft : C C[x,y], 

Q = (Qi, Q 2 , • • •, Qm) L{x) = Qi{x)Q2{x) • • • Qm{x)P{x)’' i-A- det(L(x) - y). 

Consequently, for Q = {Qi,... ,Qm) G A4 we have an affine plane curve in C^, 

given by the zeros of ifiQ). We call this curve the spectral curve. 

Each term of '^(Q) corresponds to the weight of specific highway paths as follows. We 
say that a pair of paths is noncrossing if no edge is used twice, and that a family of paths 
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is noncrossing if every pair of paths is noncrossing. Suppose p = {pi,P 2 , • ■ ■ ,Pn} is an 
unordered noncrossing family of n paths in G' using all the sources and all the sinks. The 
non-abrupt paths in p induce a bijection of a subset S C [n] with itself. We let sign(p) 
denote the sign of this permutation. The following theorem is a reformulation |MIT] . In 
our language, the proof is very similar to [LP12a[ Theorem 3.5]. 

Theorem 3.2. We have 

/(x, y) = det(L(x) - y) = ^ sign(p)wt(pi)wt(p 2 ) • ■ • wt(p„), 

P={Pl,P2,-,Pn} 

where the summation is over noncrossing (unordered) families of n paths in G' using all 
the sources and all the sinks. In other words, the coefficient x°‘y^ in f{x,y) = det(L(x) — 
y) counts (with weights) families of n paths that 

• do not cross each other; 

• cross the x-line exactly a times; 

• contain exactly b abrupt paths and n — b non-abrupt paths. 

The overall resulting sign of the monomial x^y^ is (— 

For f{x,y) = j ^ C[x,|/], we write N{f) C for the Newton polygon of 

/. This is dehned to be the convex hull of the points {{i,j) \ Oij ^ 0}. It is important 
for us to identify the lower hull and upper hull of N{f). The former (resp. latter) is the 
set of edges of N{f) such that the points directly below (resp. above) these edges do not 
belong to N{f). We exclude vertical or horizontal edges from the lower and upper hull. 

Proposition 3.3. For generic Q = {Qi ,..., Qm) G At, the Newton polygon is 

the triangle with vertices (0,n), (/c, 0) and (/c + m, 0). In iV('0(Q)), the lower hull (resp. 
upper hull) consists of one edge with vertices {k,0) and (0,n) (resp. (/c + m,0) to {0,n)). 

See § 14.11 for the proof. 


3.3. Special points on the spectral curve. For f{x,y) G C[x,y] an irreducible poly¬ 
nomial, let Gj = {{x,y) I f{x,y) = 0} C be the corresponding plane curve, and let 
G) C P^(C) denote its closure. Let Gf be the normalization of G), with a map Gf —>■ G) 
that is a resolution of singularities. We declare points on Cy- to be special if either (1) 
either x or y is 0, or (2) the point does not lie over G) (that is, x oi y is oo). 

For f{x,y) G '0(A4), dehne a polynomial gf{x,y) (resp. hf{x,y)) by f{x,y) = 
gf{x,y) other terms (resp. f{x,y) = hf{x,y) -|-other terms), where gf{x,y) (resp. 
hf{x,y)) consists of those monomials lying on the lower hull (resp. the upper hull) of 
N{f). For this f{x,y) we also dehne 

fc.= Yi f-J’ P.5) 

ij,k&Lc 


where 

Lc = {{j,i) I {m + k)i + nj 
We dehne G 0{M) by 


n(m k) — M}. 


m n'—l 

(?)=n n ^ix+jN, 

i=l j=0 


a, 


(3.6) 
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for 1 < r < iV, and define G 0{/A) by 


n 


^r{q) ■■= Ylqrj 
j=l 


(3.7) 


for 1 < r < m. 


Remark 3.4. Despite the seeming dissimilarity, (Tr{q) and er{q) have the same nature. 
Indeed, visualize variables qij as associated to crossings of two families of parallel wires 
on a torus, as it was done in Section [2j Then ar{q) is the product of parameters on 
the r-th horizontal wire (out of N), while er{q) is the product of parameters on the r-th 
vertical wire (out of m). In particular, the symmetry between Q and Q from Section [2] 
switches the ar{q) and the er{q) into each other. 

Lemma 3.5. For f{x,y) G we have 

N 


9t(x, y) = + ffr A), 


(3.8) 


r=l 




(3.9) 


See § 14.21 for the proof. 

In the following, for f{x,y) = '0(<?), q E Ai we study the special points of Cf related 
to polynomials gf{x,y), hf{x,y) or /(x, 0). 

We have f{x, 0) = det P{xY Y\T=i Qi{x), where det Qi{x) = 11^=1 + (—= 

ej(g) + (—Thus the roots of /(x,0) = 0 are exactly the (—l)”ej,(g), where er{q) 
is dehned by fl3.7p . It is also clear that ei,..., depend only on f{x, y). The following 
lemma is obtained immediately. 

Lemma 3.6. Suppose that ei,..., are distinct and nonzero. Then there are m special 
points Ai = ((—l)”ei(g), 0) on Cj with y = 0 and x is nonzero. 

The point (0, 0) on C'f is usually singular, whenever k >2. 

Lemma 3.7. Suppose cxi,..., ctat are distinct and nonzero. Then there are N points of 
Cf lying over (0,0) G Cj. 

Proof. By fl3.8lh the meromorphic function y"^'/x^' takes the N distinct values 
for r = 1, 2,..., as (x, y) —)■ (0, 0). So there are at least N points on Cf. On the other 
hand, looking at N{f) we see that analytically near (0,0), the polynomial / can factor 
into at most N pieces. □ 

Let Oi, O 2 , • • •, Oat denote the the special points of Lemma 13.71 Near Or there is a 
local coordinate u such that 



It turns out that Cf C P^(C) has only one point at 00 . Due to the polynomial hf{x, y), 
in homogeneous coordinates, this point is 


[1:0:0] a n > k + m, 


P' = < [1 : 1 : 0] a n = k + m, 
I [0 : 1 : 0] if n < fc + m. 
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To compute this, we first homogenize f{x,y) to get F{x,y,z) = f{x/z,y /z)'wheTe 
d := max(n, m + k). Then we solve F{x, y, 0) = 0. Recall that fc was dehned in fl3.5p . 

Lemma 3.8. Suppose that fc 7 ^ 0. Then there is a unique point P ^ Cf lying over P'. 

See § 14.31 for the proof. 

3.4. A good condition for the spectral cnrve. Let Vn^m,k be the subspace of C[x,y] 
given by 

n—1 


Vn,m,k = I ((-!/)" +a: + ^y^fi{x) 


i=0 


M^) = e C[a;] J> , (3.10) 

71,771, A; 

where we write Ln,m,k for the set of lattice points in the convex hull of {{k, 0 ), ( 0 , n), {k + 
m,0)} but not on the upper hull. By Proposition 13.31 and Lemma 13.51 we have that 
(Z Vn,m,k- 

Definition 3.9. Define the subset V C Vn^m,k as the set of f{x,y) G Vn^m,k, satisfying 
the conditions: 

(1) f{x, y) is irreducible; 

( 2 ) C'f is smooth; 

(3) fc ^ 0; 

and such that the special points on Cj consist exactly of: 

(4) m distinct points Ai, A 2 ,..., Am where := ((—l)”ej, 0) and e* 7 ^ 0; 

(5) N distinct points Oi, O 2 , ■ ■ ■, On lying over (0, 0), where near Or there is a local 
coordinate u such that 


and (Jr 7 ^ 0 ; 

( 6 ) a single point P lying over the line at infinity P^(C) \ C^. 
For / G V, the genus g of Of is given by 


(3.11) 


9 = ^ ((n- l)m 


M-A + 2). (3.12) 

Indeed, it follows from Pick’s formula and Proposition 13.31 that the number of interior 
lattice points of N{f) is equal to the right hand side. This formula for the genus then 
follows from |Khol Corollary on p. 6 ]. Alternatively, the genus can be computed using 
the Riemann-Hurwitz formula, as in |vMM] . 

Proposition 3.10. For f ^ V, we have that f3~^{f) 7 ^ 0. Moreover, the set V := 
Vn'^(A4) is a Zariski-dense subset of'ip^Ai). 

We give the proof in § 14.41 Informally, the last statement of Proposition 13.101 states 
that most curves in satisfy the “niceness” conditions listed in Definition 13.91 

From Definition 13.91 it follows that for / G V, the meromorphic functions x and y on 
Of satisfy 

N N m 

{x) = n''^Oi - nP, (y) = A' ^O* + ^ A* - (m + k)P. 


2=1 


2=1 


2=1 


(3.13) 
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Remark 3.11. The condition that Ai, ..., (resp. Oi, O 2 , ..., Ojv) are distinct imply 
that the qnantities ei,..., (resp. ai,..., a^) are distinct. Many of onr main resnlts 
still apply after a modihcation even when this condition does not hold. 


4. Proofs from Sections [2] and [3] 

4.1. Proof of Proposition [H731 We use the interpretation of /(x, y) given by Theorem 
13.21 Let p = (pi,... ,pn) be a family of noncrossing highway paths in G', as in Theorem 
and let wt(p) = wt(pi) • ■ •wt(p„). 

There are k + m opportunities for a highway path in p to pick up the weight x: from 
each of the m vertical wires and from the k horizontal wires crossing x-line. For a hxed 
power (where 6 = 0,1,..., n), or equivalently a hxed number of abrupt paths, we will 
bound the maximal and the minimal possible value of the exponent a of x. 
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Figure 6. The red path crosses x-line less than the green one. 

The hrst key observation is that we can consider p to be a family of noncrossing closed 
cycles on the toric network G. An abrupt path p is simply the “cycle” that starts and 
ends at the vertex i (= source i and sink i identihed) in G and does not move. 

Lift such a highway cycle G to the universal cover, as shown in Figure [6l We obtain a 
path that starts and ends at two vertices labeled with the same integer i G {1, 2,..., n}, 
but £ periods (of m vertical lines each) to the right of the original source. More precisely, 
if G is obtained by gluing paths Pij,... ,p^ in G', then G has length £. The x-line now 
has a staircase-like shape as shown in Figure |6l We claim that if G has length i, then it 
crosses the x-line at least £/c/n times and at most i{k + m)/n times. 

To see this, observe that if we lift the ending vertex several periods (consisting of n 
horizontal wires) up, we preserve the length and we increase the number of crossings with 
the X-line. Similarly, if we lower the ending vertex several periods down, we preserve the 
length and we decrease the number of crossings of x-line. Thus, the smallest and the 
largest ratios of the quantity (crosses of x-line/length) is achieved for the lowest and the 
highest possible highway paths, shown in Figure [6] in purple and green. The former is 
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the horizontal path, while the latter is an alternating right-np staircase path. For these 
paths, the ratios are exactly k/n and {k + m)/n. 

4.2. Proof of Lemma 13.5L One of the conseqnences of the proof of Proposition 13.31 is 

that the monomials for which the lower bonnd k/n oi the ratio is reached are the ones 
coming from horizontal highway paths on the nniversal cover. Let ns call each snch closed 
cycle a coil. For example, the pnrple line in Fignre E] represents a coil passing throngh 
sonrce i, as well as passing throngh the n' other vertices between 1 and n that have 
residne i modnlo N = gcd(/c,n). Thus the terms of g{x,y) are formed in the following 
way: for each of the N coils we decide whether to include it into our family of paths, or 
to make all paths starting at its sources abrupt. The second choice corresponds to the 
contribution The hrst choice gives (Hlii 11^=1 Qi,r+jN)x’^\ which is the weight of 

that coil. Thus, the r-th coil contributes the factor of ((nr=i n;li qi,r+jN)x^' + (-?/)”'), 
and fl3.8p follows. 

Similarly, the monomials for which the upper bound {k + m)/n of the ratio is reached 
are the ones coming from right-up staircase paths on the universal cover. Let us call each 
such closed cycle a snake path. For example, the green line in Figure [6] represents a snake 
path passing through source i, as well as through all the n/M other vertices between 1 
and n that have residue i modulo M = gcd(m + k,n). Thus the part contributing to the 
upper hull of f{x, y) is a product of factors ^ where the term {—y)'^^^ 

corresponds to choosing to have abrupt paths, while corresponds to choosing 

to have the snake path. Thus we obtain fl3.9p . 

4.3. Proof of Lemma 13.81 We analyze the singularity at P' G P^(C) on the chart 
1 / 7 ^ 0, that is, {[x : 1 : z] \ x^z E C}. Let Cj be the affine curve given by f{x,z) : = 
F{x, 1, z) = 0. 

(i) n < k + m: using 03.91) we can write f{x, z) as 

max(m,m+A:—n) 

f{x,z) = — ^ Z^fj{x) 

i=i 

where deg^fj{x) > 1. Then, when m + k — n = 1, is smooth at (0,0) since 
df /dz\{Q^o) = (—1)”/^ 7 ^ 0. When m + k — n > 2, the point (0, 0) is singular. Dehne K := 
C{x)[z]/{f{x, z)). We consider all valuations n : iF -» Z on iF satisfying that v{x) > 0 
and v{z) > 0. We shall show that such valuations are unique, from which the uniqueness 
of P follows. Dehne hf{x,z) := {{—1 )^z^^'m " which is the part of f{x,z) 

consisting of monomials lying on the lower hull of N{f), corresponding to the upper hull 
of 7V(/). Then we see th.a.t v{f{x, z)) = v{hf{x, z)) = M ■ v{x^^ + {—1)^z^^ ). To 
be consistent with v{f{x,z)) = n(0) = oo, the only choice we have is v{x) = ^ and 
v{z) = ^. 

(ii) n = k + m: we can write f{x, z) as 

m 

k=l 
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where /^(x — 1) G C[x — 1]. Thus P' is smooth point on Cj if /i(0) 7 ^ 0. Looking at the 
Newton polygon N{f) and ([33]), it follows that /^(O) = = fc ^ 0- 

(hi) n > k + m\ this case is nearly the same as the case n < k + m. 

4.4. Proof of Proposition [37TOl The hrst statement follows from Theorem 15.31 in the 
next section. In the following, we prove the second statement. 

We hrst show that V contains a Zariski-dense and open subset of A Zariski- 

open subset of consists of /(x, y) with Newton polygon iV(/) given by Proposition 

13.31 This polygon is not a non-trivial Minkowski sum of two other polygons, so f{x,y) 
is irreducible. Similarly the conditions that C'f is smooth and fc^O are Zariski-open 
conditions on Vn,m,k- By Lemma 13.81 fc ^ 0 implies that ( 6 ) in Dehnition 13.91 holds. 
The calculations in the proofs of Lemmas 13.61 and 13.71 then imply that V contains a 
Zariski-open and dense subset of 14,m,fe¬ 
lt thus suffices to show that contains a Zariski-dense subset of Vn,m,k- We shall 

use the following result. 

Lemma 4.1. The set a{A4) is Zariski-dense in C. 

Proof. This follows from [LPlll Proof of Theorem 4.1] where it is shown that the map 
a \ M. ^ C is, generically a m! to 1 map between two spaces of dimension mn. □ 

By the hrst statement of Proposition [3301 we have V C /3(4). By Lemma 1431 4(34) = 
/3(a(M)) contains a Zariski-dense subset of V. It follows that 4(34) is Zariski-dense in 
Vn^rn,k- This Completes the proof of the second statement of Proposition 13.101 

Remark 4.2. The nonconstant coefficients /jj of f{x,y) can be pulled back to func¬ 
tions on 34 or £. A consequence of our proof is that 4(34) is Zariski-dense in I4,m,fe, 
and therefore the functions fij are algebraically independent on 34 or C. It would be 
interesting to obtain a direct proof of this. 

4.5. Proof of Theorem 12.41 We employ the technique hrst introduced in |LP13] . 



Figure 7. Crossing merging and crossing removal moves with vertex 
weights shown. 




Figure 8 . Yang-Baxter move with transformation of vertex weights shown. 
Specihcally, we realize the geometric i?-matrix transformation, dubbed the whirl move 
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in |LP13j . as a sequence of local transformations on our toric network. The trans¬ 
formations we shall employ are of three kinds: crossing merging/unmerging, crossing 
creation/removal, shown in Figure O and Yang-Baxter move shown in Figure [HI 

The whirl move R occurs between two parallel wires adjacent to each other and wrap¬ 
ping around a local part of the surface that is a cylinder. The way it is realized as a 
sequence of local moves is illustrated in Figure [HI First a crossing is created with weight 
0, and split into two crossings of weight p and —p. One of them is pushed through the 
wires crossing our two distinguished wires, until it comes out on the other side. As it is 
proven in |LP131 Theorem 6.2], there is at most one non-zero value of p for which the 
end result is again a pair of crossings of weights p and —p, and thus those two can be 
canceled out. The resulting action on the weights along two parallel wires is exactly the 
whirl move i?, and does not depend on where the original auxiliary crossing was created. 
The value of p does depend on the location j where the new crossings are created, and 


XX 


yi 
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y 2 

^3 


ys 








Xi 


yi 

X 2 


y 2 

X3 


ys 







Figure 9. For a unique choice of the weight p, the weight that comes out 
on the other side after passing through all horizontal wires is also p; the 
resulting transformation of x and y is exactly the whirl move. 


is given by 


Here the variables Xi and Pi are weights along the two wires as in FigurelHl and 

is the energy of the cyclically shifted vectors x^^^ = {xj+i, Xj+ 2 , ■ ■ ■, Xj) and y^^'> = 


{Pj+i, yj+ 2 , • • •, Vj) as introduced in § O 


Now, the fact that R satishes the braid move, that is, and 

S£Si+iSi = S£+iSiS£+i is shown in |LP131 Theorem 6.6]. Indeed, these relations happen 
at a local part of the surface (in our case, torus) that looks like a cylinder. 

On the other hand, the commutativity of the si and the S£ does not follow from |LP131 
Theorem 12.2]. This is because in |LP13j we only considered the case when the pairs 
of parallel wires intersect once. In our case on the torus however, it is common to have 
horizontal and vertical wires intersect more than once: this happens for any h 7^ 1. The 
proof in such a situation is essentially the same. 

Indeed, if we have two pairs of parallel wires crossing as in Figure UHl but possibly 
more than once, we can realize each of the two corresponding i?-moves by a sequence 
of local moves as above. It is a local check that performing one of the two sequences 
does not change the value of p = one needs to perform the other, because each 

of the n Xi, n Pi and E{x^^\y^^'^) is not changed. It is also a local check that the two 
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Figure 10. 

sequences commute once the parameters p for each are chosen, this is [LP131 Proposition 
3.4]. Thus, commutativity follows. 

5. Eigenvector map 

In this section, we hx / G V (see Dehnition l3.9p and consider the corresponding smooth 
curve Cf. For j G Z, we set Oj := if j = r mod N. Similarly, for j G Z, we set 
Aj := Ar ii j = r mod m. 

5.1. Generalities. A divisor D = on an algebraic curve G is a hnite formal 

integer linear combination of points Pi on C. We write H > 0 if > 0, and say that D 
is positive in this case. The degree of D is given by deg(Zl) = 

Given h a meromorphic function on C, we let {h) = {h)o — {h)oo be the divisor of 
h. Here, (h)o denotes the divisor of zeroes, and (h)oo denotes the divisor of poles. Two 
divisors Di, D 2 are linearly equivalent, Di ~ ZI2, if there exists a meromorphic function 
h such that {h) = Di — D 2 . We write [D] for the equivalence class of D with respect 
to linear equivalence. The Picard group Pic(G) is the abelian group of divisors on C 
modulo linear equivalence. For j G Z, we write Pic-^(G) for the part of the Picard group 
of C that has degree j, that is, Pic^(G) ;= {D : a divisor on C \ deg{D) = j}/ ~. 

To each divisor D we associate a space of meromorphic functions 

C{D) = {f\{f) + D>0}. 

If D = — 0; then in words, C.{D) consists of meromorphic functions which are 

allowed to have poles of order rii at Pi. We have dim(£(Zli)) = dim(£(Zl2)) if Di and 
D 2 are linearly equivalent. 

5.2. Positive divisors of a Lax matrix. A divisor D G Pic®(G/) of degree g is called 
general if dim£(Zl) = 1. A divisor D G Pic(G/) is called regular with respect to the 
points P and Oj if D is general and if dim£(Zi) + kP — J2]=n-k^j) = 0 for fc > 0. 
For brevity, we will sometimes just say that a divisor is regular when it is regular with 
respect to P and the Oj. 

Let us now hx L{x) G /3“^(/). Let 

A.J := (-l)‘+'|L(l) - yUj 
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denote the (signed) (i, j)-th minor of the matrix L{x) — y. 


Theorem 5.1 (cf. |vMM] ) . There exists a positive divisor R of degree 2g+2 supported on 
Sp,o ■= {P, Oi) • • •) On}, and uniguely defined positive general divisors Hi, H2,..., -D„, Hi, 
of degree g such that for all {i,j) G [n]^, we have 

i—1 

i^hj) ~ Oj + Di + {j — i — 1)P + Or — R- 

r=j+l 

In addition, Hi, ..., H„ have pairwise no common points, and Hi, ..., H„ have pairwise 
no common points. 


See § 16.41 for the proof. Note that Yll=j Or is to be interpreted in a signed way: 


i 00 

:= 

r=j r=j 




( Y!r=j Or i > j, 

< 0 * = J -1, 


For L{x) G fi~^{f), dehne gi ;= Since L{x) — y is singnlar along Cf, the vector 
g = {gi,g 2 ,... ,gn)~^ (thonght of as a vector with entries that are rational fnnctions on 
Cf) is an eigenvector of L{x) — y. We dehne g^ for i G Z, by gi+n = x~^gi. We also dehne 
K '■= gi/gn for i G Z. Thns = 1 and hi+n = x~^hi. The vector (hi, ^ 2 ,.. ., is 
also an eigenvector of L{x) — y. 


Definition 5.2. For L{x) G fi ^(/), dehne the divisor V = T>{L{x)) on Cf to be the 
minimnm positive divisor satisfying 

n 

{hi) +V> 

j=i+l 

that is, hi G C{V + {n — i)P — E^=i+i Oj), for i = 1,..., n — 1. 

It follows from Theorem 15.11 that H = H„ from the theorem and that this divisor is 
nniqnely determined by Dehnition 15.21 In particnlar, H is a positive regnlar divisor of 
degree g with respect to the points P and Oj, and we have 

n 

(fo) = D,-V-{n-i)P+Y, Oy (5.1) 

j=i+l 

5.3. The eigenvector map. Let Ra (resp. Rq) denote the set of orderings of the m 
points Ai,...,Ara (resp. Oi,..., On)- 

Ra := {^{Ai,..., Am) I i' G 6m}, Ro ■= {h{Oi,..., On) \ h G &n}- 

By onr assnmption that / G V, the points Ar (resp. Or) are distinct, so Ra has cardi¬ 
nality ml and Rq has cardinality N\. 

We write Ur {r = 1,... ,m — 1) (resp. fo (r = 1,..., iV — 1)) for the generator of 6m 
(resp. 6iv), which permntes the r-th and the (r-|- l)-st entry of A G Ra (resp. O G Rq)- 
By onr assnmption that / G V, we have fcfi^ 0- Dehne Sf C (C*)^ by 
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For i = 1,..., M, define q to be the coefficient of the maximal power of x in the 
{i,i + 1) entry of . (Here, maximal power is the maximal power for a generic 

L{x) G C.) An explicit formnla of q for L{x) G a{A4) is given in Lemma [6.51 
We now dehne the eigenvector map pf for / G V by 

Vf ■■ r'(/) ^ WiCf) xSjxRo (5.2) 

L(x) I —> {[V],{ci,...,cm),0), (5.3) 

where the ordering O = {Oi,... ,On) is uniqnely determined from (Ajj) by Theorem 
15.11 For / G V we dehne the map </>/ by 

0/ : 0"'(/) ^ W{Cf) xSjxRoxRa (5.4) 

Q = (Qi, (ci,..., Cm), O, A), (5.5) 

where A = (((-l)”ei, 0),..., ((-l)”em, 0)), and er{q) is given by fl^ . That (ci,..., cm) 
lies in iSj is the content of Lemma [6.61 below. We will omit the snbscripts of and (pf 
and just write rj and 0 when no ambiguity will arise. 

The following theorem should be compared to Theorem 1 in [vMMj . 

Theorem 5.3. We have a one-to-one correspondence between L{x) G /3“^(V) and the 
following data: 

(a) / G V, 

(b) (["D], c, O) G Pic^(C/) X Sf X Rq where V is a positive divisor on Cf of degree g, 
and regular with respect to P and the Oj. 

See 1 16.51 for the proof. It is shown in |LP11( Proof of Theorem 4.1] that the map 
a : M —)■ £ is generically m! to 1. Let Xi denote the open subset of Ai where (a) the 
map a is ml to 1 (that is, Ai = a{M.) is a m! to 1 map), and (b) the 

hnite symmetric subgroup acting via the i?-matrix is well-dehned on AA. 

Theorem 5.4. Fix / G V. We have an injection from fj~^{f) A Ai to the collection of 
data ([P],c, 0,A) G Pic^(C'/) x Sf x Rq x Ra such that 

(a) V is a positive divisor on Cf of degree g. R is regular with respect to P and Oj-s. 

(b) c = (ci, ...,cm)^ Sf. 

(c) O G Ra¬ 
id) A G Ra- 

See 1 16.61 for the proof. We will usually just write fj~^if) instead of 0“^(/) AM. when 
no confusion arises (for example, when discussing the action of the i?-matrix on the 
spectral data). 

Let us perform a quick dimension count. The dimension of AT or £ is equal to 
mn. The dimension of Pic®(C'/) x 5/ is equal to 5 ^ + M — 1. The number of lattice 
points in the convex full of {{k, 0), (0, n), (/c + m, 0)} is equal to g + N -\- m -\- M. Thus 
the dimension of Vn^m,k f|3.10p is equal to g + N + m — 1. Using 03.121) we obtain 
dim(AT) = dim(Pic^(U/) x Sf) + dim(14,m,fe), consistent with Theorems 15.31 and 15.41 
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6. Proofs from Section [5] 

The first three subsections are devoted to introduce key lemmas for Theorem 15.11 and 

[531 

6.1. Special zeroes and poles of the eigenvector. 

Proposition 6.1. The eigenvector g = {gi)f of L{x) G satisfies the following. 

(i) For any j G Z, the rational function gj/gjj^i on Cf has a zero of order one at 

Oj+i. 

(ii) For any j G Z, the rational function gj/gj+i on Cf has a pole of order one at P. 

To prove this proposition, we shall need the following generalization of Theorem 13.21 
whose proof is essentially the same as that of Theorem 13. 2 1 Suppose p = {pi,P 2 , ■ ■ ■, Vn-i} 
is an unordered noncrossing family of n paths in C using all the sources except i, and all 
the sinks except j. Identifying [n] \ i ~ [n — 1] ~ [n] \j via the order-preserving bijection, 
we have that the non-abrupt paths in p induce a bijection of a subset S' C [n — 1] with 
itself. We let sign(p) denote the sign of this permutation. 

Let \{L[x) — y)\ij denote the maximal minor of {L{x) — y) complementary to the 
(i, j)-th entry. 

Theorem 6.2. We have 

\{L{x) - y)\ij = sign(p)wt(pi)wt(p 2 ) • ■ • wt(p„_i), 

P={Pl,P2,-,Pn-l} 

where the summation is over noncrossing (unordered) families of n — 1 paths in C using 
all the sources except i, and all the sinks except j. In other words, the coefficient of x°‘y^ 
in |(T(x) — y)\i,j counts (with weights) families of highway paths that 

• start at all sources but i and end at all sinks but j; 

• do not cross each other; 

• cross the x-line exactly a times; 

• contain exactly b abrupt paths and n — b — 1 non-abrupt paths. 

Example 6.3. Let (n, m, k) = (3, 3,1). The Lax matrix is given as follows: L{x) = 

( (hp(l2,W + gi, 1 ^ 3 , 33 ; + ^ 2 , 3 ^ 3 , 33 ; ^ 1 , 1 ^ 2 , 1 ^ 3 , 13 ; X^ 

51,2^2,2^3,2 + X qi,2X + q2,lX + qi,2q2,2X + qi,2q3,lX + q2,iq3,lX 

V 91,3^2,3 + 91,3^3,2 + ^2,253,2 91,3^2,3^3,3 + X gi^gX q2,2X + q^pX 

Consider the minor \{L{x) — |/)|i ,2 = with rows 2,3 and columns 1,3. Here are some 
terms that appear in it: 

|(L(x) - I/) 1 1,2 = q 2 , 2 x‘^ - qi,3q2,3qi,2q2,2X -xy - . 

The term q 2 , 2 X^ for example is formed by two paths: one starting at source 3, turning 
at gi, 3 , turning at qip, going through at ^ 2 , 2 , turning at turning at q^p and thus 
hnishing in sink 3; the other a staircase path starting at source 2, turning at each of gi, 2 , 
gi,i, g 2 ,i, q 2 , 3 , q 3 , 3 , q 3,2 and hnishing at sink 1. The hrst path contributes weight q 2 , 2 X, 
while the second contributes weight x. Note that source 1 and sink 2 remain unused, 
as they should according to the theorem. Also note the the paths give bijection 2 1 —>■ 1 
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and 3 I—)■ 3 between the used sources and sinks. The induced permutation is the identity 
permutation, and that is why we have sign + in front of this term. 

The term —xy corresponds to one abrupt path of weight —y from source 3 to sink 3, 
and one staircase path, the same as for the previous term. The term —gi, 3 ^ 2 ,sQ'i, 2 ^ 2 , 23 ^ 
corresponds to two paths that induce the map 2 1 —)■ 3 and 3 1 —)■ 1 on sources and sinks, 
of weights qifiq 2 , 2 X and qi^ 3 q 2,3 respectively. The minus sign arises since the induced 
permutation in S 2 is a transposition. 

Let us prove Proposition 16.11 (i) If = 1, the result is easy: gj/gj+n vanishes to order 
n at Oi- We also have a shift automorphism s : C ^ C (see the proof of Proposition 
Ib.lf iiP which sends Oi to Oi and pulls gj/gj+i back to gj+i/gj+ 2 - We will thus assume 
iV > 1. 

For each i, dehne 

yW ;= ((-l)i|L(x) - y\i^i ,..., (-l)"|L(x) - y\i,n)^ . 

We shall think of as a vector whose entries lie in the coordinate ring of the affine plane 
curve Cf. Like the vector g = {gi, g 2 ,, gn)^, the vectors are (nonzero) eigenvectors 
of the matrix L{x). The matrix L{x) is singular along the curve C*/, but generically it 
has rank n — 1. Thus for any i, the vectors g and are multiples of each other. More 
precisely, g/v^'^'^ is a rational function on Cf. To show that gj/gj^i has a zero at 
we shall calculate using a convenient choice of 

Choose V = Then Vj = \L{x) — y\j^j and Vj+i = \L{x) — Thus, Vj counts 

the families of paths that start at all sources but j and end at all sinks but j, as in 
Theorem 16.21 

Let us make the substitution (x, y) u^') inside Vj, and let Vf{x,y) 

denote the terms of in Vj{x,y) that give the lowest degree in u after the substitution. 
Call this degree d. The terms in Vj{x,y) are obtained by either taking abrupt paths, or 
by taking coils. In the proof of Lemma 13.51 we defined the N coils in the network G, 
which we denote Ci, C2, ..., Cat. To obtain a path family p contributing to n', instead of 
Cj, we include the n' — 1 abrupt paths which use the vertices j' where j 7 ^ j' but j = j' 
mod N. For each of the coils Ct where t ^ j, we can either include the coil (that is, 
include the n' paths in C) in p, or we use the corresponding n' abrupt paths instead. In 
particular, considering Cj+i we see that n' has a factor of [aj+ix^' + (—?/)"''). This factor 
vanishes under the substitution {x,y) ~ {u^\ — u^'), thus creating a zero at 

the point Oj+i of order at least d + 1. 

For Vj+i, we count families of paths that start at all sources but j and end at all sinks 
but j + 1. Let {x,y) denote the terms of lowest degree in Vj+i{x,y). This lowest 
degree is again equal to d as well (we caution that if j = n, then we take Vj+i := x~^vi). 
The calculation of Vj_^_i{x,y) is similar to that of Vj{x,y), except that instead of a single 
“incomplete” coil with index j modulo N, we have two incomplete coils Cj and Cj+i 
with indices j and j + 1 modulo N . We obtain 

v'j+i{x,y) = a(g)a:“/ JJ {atx^' + (-?/)”') 
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where a{q) is some nonzero polynomial in Qab-s, and a,f3 are nonnegative integers, and 
the prodnct is over t G {1, 2,.. . N} not eqnal to j or j + 1 modnlo N. As a resnlt, we 
see that Vj+i vanishes at Oj+i of order exactly d. 

Thus Qj/gj^i has a zero at Oj+i- The fact that gj/gj+i vanishes to exactly order 1 
follows since we know that gj/gj+n vanishes to the order of n' at each Oj. 

(ii) For L{x) G /9“^(/), we have rational functions gi/gi^i on Cf. To emphasize the 
dependence of gi/gij^i on L, we write gijgiJ^i^L). Let be the shift map given 

by L{x) i-A- P{x)L{x)P{x)~^. (We will introduce the shift operator ^ on AT in 11 17.21 1 
Then (]{(;*{L{x))) = P{L{x)) is invariant under r. Let F be a Zariski dense subset of 
P~^{V) such that = Y. Then there exists a nonnegative integer a* given by 

tti = max{oTdp{gi/gi+i{L))^ | L G F}, 

where OTdp{gi/ gipi{L))^ denotes the order of the pole of gi/gipi{L) at P. Since gi/gipi{g*{L)) 
gi_i/gi{L) as rational function on we conclude that for all i, we have ordp^gi/gipi{L))ao 

a := minj{aj}. But we know that ordp^gi/ gipn{L))oo = R for all i and L E C. Thus we 
must have a = 1. Since this value does not depend on the choice of F, the claim follows. 

6.2. A holomorphic differential. Let C, be the differential form on the curve Cf given 
by 

^ x^~^dy 
^ = 9 / ■ 

dx 

Lemma 6.4. The divisor of the differential form ( is supported on Sp^o = {P) Oi ,..., Ov}-' 

N 

(c)= {k' -i)y: Oi + ((n — l)m — k — M) P. 

i=l 

When g >1, it is holomorphic. 

Proof. By using local expansions of / around Oi and P, we get 

= {n' — k'n) ^ Oi + n{k + m — 1)P, 

' i=l 

N 

{dy) = {k' - 1) E Oi — {m + k + M)P. 

i=l 

From these and {x) in §4.3, we obtain {(). 

Now we prove that if > 1 then {n — l)m — k — M > 0. Let p be the number of 
integer points inside a parallelogram in M^, whose vertices are (0,n), {k,0), {k + m, 0) 
and (m, n). We have p = 2g + M — 1, since the parallelogram is composed of the Newton 
polygon N{f) and its copy, sharing the upper hull of N{f). Then the claim is equivalent 
to that if > 1, then p + N — l>M + k. 

When k = n, N = n follows. Then we have p + N — 1 — (M + k) = p — {M + 1). It 
reduces to 2g — 2, which is non-negative when g > 1. 

When m + k = n, M = n follows. Then we have g = |(n(m — 1) — m — A -|- 1) 
which is negative if m = 1. So we assume m > 2. We have p + N — 1 = m{n — 1) and 
M + k = 2n — m. Thus, we obtain p + N — 1 — (M -|- fc) > 0 when m > 2. 
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When n > k, we prove the claim by choosing M + k different points from p points in 
the parallelogram. When m + k < n (resp. m + k > n), we can choose M — 1 points on 
the npper hnll of N{f), k points; (i, n — i) for i = 1,..., A; inside the npper (resp. lower) 
triangle, and one point inside the lower (resp. npper) triangle. 

Finally the claim follows. □ 

6.3. About (ci,..., Cm) G Sf. Recall that q is the coefficient of the maximal power of 
X in the (£, £ + 1) entry of dehned in 1 15.31 

Lemma 6.5. When L{x) = a{q = (g^)), we have 

ce = Qij- 

i+j=i+l mod M 

Proof. By Lemma l4.ll it snffices to prove the lemma for f{x,y) G We apply 

Lemma [3. II to interpret the q as connting almost snake paths in G, that is, paths that 
tnrn at all steps bnt one. The almost snake paths p that are connted start at vertex i and 
end at vertex £+1, and when drawn in the cylindric network G', they are disjoint nnions 
of n/M paths. In other words, p wraps aronnd the pictnre Fignre 0] horizontally n/M 
times. Snch paths p are completely determined by the single vertex that the path goes 
throngh. With i hxed, the weights of these vertices are exactly the qij with i+j = 1 + 1 
mod M. □ 

Recall that we dehned /c in fl3.5l) . 

Lemma 6.6. When L{x) G we have 

M 

JJq = /c = ± Ai- (6-1) 

£=1 {j,i)eLc 

Thus the product of all ci is constant on 'ip~^{f). 

Proof. We hrst show that 

degree of fij in the qr,s = (m + k)n — ((m + k)i + nj). (6.2) 

Indeed, lift a family p of paths that contribntes to fij (as in Theorem l3.2l) to the nniversal 
cover, as in Fignre [6l By Theorem 13.21 the family p inclndes exactly i abrnpt paths and 
n — i non-abrnpt paths. Let ns snppose that the sonrces used by the n — i non-abrupt 
paths are (1/2, zi),..., (1/2, Zn-i) and the sinks are (m + 1/2, tci),..., (m + 1/2, Wn-f), 
where 1 < Zr < n and 1 < Wr < n + m. To agree with our convention in Figure 0] that 
source and sink labels increase as we go down, we will take the y-coordinate to increase 
as we go down in Figure [HI but otherwise, the coordinates we are using are the usual 
Cartesian coordinates. 

For each r, let tcj. G [1 — k,n — k] be chosen so that w'^. = Wr mod n. From Theorem 
13.21 it also follows that 

n—i 

^(tCr — wl)/n = i, 

r=l 

as this is the number of times the x-line would be crossed by the paths. Also, we 
know that {tci,..., Wn-i} = {+ — k,..., Zn-i — k} mod n and thus {tcj,..., = 
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{zi — k,..., Zn-i — k}. This is because our collection of non-abrupt paths must form a 
collection of cycles in G. The number of the qr,s a path from {1/2, Zr) to (m + 1/2, Wg) 
picks up is equal to Zr + rn — Wg, and summing over all the non-abrupt paths in p we get 


n—i n—i 

degree of fij in the gj.,s = m(^n — j) + Zr — Wg 

r=l s=l 

n—i n—i 

= m{n - j) + ^Zr - ^w',+ jn 

r=l s=l 

= m{n — i) + {n — i)k — jn = {m + k)n — {{m + k)i + nj). 


as claimed. This degree is equal to M when (j, i) G L^,. 

Now we prove the lemma. Let us partition all the edges in the network G into M 
snake paths, Hi,..., Hm- Suppose p is a closed path in G. Then we can partition p 
into snake intervals - maximal contiguous segments of p that follow one of the H^. Such 
snake intervals are connected by a horizontal step going through a vertex (and picking 
up its weight), and moving from a snake path Hi to the next snake path It is easy 

to see that the vertices which separate Hi from iL^+i are exactly the ones labelled with 
qtj where i + j = ^ + 1 mod M. 

In order for p to be closed, it must consist of c horizontal steps, where c is a multiple of 
M, since p must come back to the snake path it started at. When (j, i) G L^,, 06.21) shows 
that all families of toric paths that contribute to fij consist of just a single closed path 
in G that picks up exactly one of the qij from each of the sets Ti = {qij \ i+ i = i+ l 
mod M}, for ^ = 1,..., M. Using Lemma 16.51 we see that each such term appears 
exactly once in the product 

We need to show that this is not only an injection, but a bijection. For each term 
contributing to G-, h suffices to hnd a closed path in G that has this weight. By 
06.2p . such a closed path would necessarily contribute to fij for (j, z) G L^,. 

For each appearing in our chosen term of G, we draw a horizontal step through 
the vertex labelled in the network G. From the endpoint of each such horizontal step 
(say from Hi to attach a snake interval in leading to the start point of the 

horizontal step which goes from iL^+i to iLf+ 2 - Some of these snake intervals may be 
empty. When we glue everything together, we get the desired closed path. □ 


6.4. Sketch proof of Theorem 15.11 The proof is analogous that of [vMMj . We explain 
the differences. When n — k < m < 2n — k, the triple of integers [N, M, M') in van 
Moerbeke and Mumford’s work [vMM] corresponds to the parameters (n, n — k, m+k—n) 
in our case. 

One hrst shows that there is a positive regular divisor of degree g, denoted V, satisfying 
Dehnition 15.21 The correspondence L{x) ha P is the main construction in |vMM[ Theo¬ 
rem 1]. Their results do not formally apply since our Lax matrix is not regular in their 
terminology. Nevertheless, the properties of L[x) ha P is proved in the same manner 
as |vMMl Lemmas 3 and 4], where Proposition 16.11 takes the place of |vMMl Lemma 2]. 
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The divisors Di and Dj of Theorem 15.11 are obtained by shifting and transposing the 
Lax matrix. [vMMi Proposition 1] then says, in our terminology, 

i—1 

(A,,.C) = +D[ + {j-i-l)P+Y. 

r=j+l 

where C, is the differential of Lemma [6.41 We now take i? = (C) to obtain Theorem 15.11 
The statement about common points is on [vMMl p.ll7]. 

6.5. Proof of Theorem 15.31 We shall also use the following lemma which is proved in 
a similar way to |vMMt Lemma 5]. 

Lemma 6.7. Suppose V is a positive divisor on Cf of degree g. which is regular with 
respect to the points P and Oj. We have 

(i) dim C{V + {n - i)P - YJj=i+i Oj) = 1 for i e Z. 

(ii) Suppose for each i E Z, we have fixed a nonzero element hi G C{V + (n — i)P — 
Y!j=i+iOj)- Then for ii < i^, and any h G >^(1) + (n - ii)P - Xl”=* 2 +i 

are unigue scalars hi^^,..., bi^ such that 

i2 

h = ^ bihi, 

i=ii 

with bi^, bi^ 7 ^ 0. 

Let us prove the theorem. Due to Theorem 15.11 Lemma 16.61 and the assumption on V, 
for L{x) G /S“^(V) we have / := l3{L{x)) G V, and rif{L{x)) = ([D],c, O) satisfying (b). 
In the following we construct the inverse of ? 7 / fl5.2p for / G V. Around P E Cf we have 
the local expansion {x, y) = (xom”"', you~"^~^) as n —)■ 0. We take c = (ci,..., cm) G Sf, 
and recursively dehne di {i E Z) by 

dn = 1, di = Q dj+i for i = f mod M. 

Let hi E C{V+ {n — i)P — Yf'j=iJ^i Oj) have an expansion around P as hi = diU~"'^^ + • ■ ■. 

Since yhiE C{V + {n + m + k — i)P — Yl^j=-k+i+i Oj)) using Lemma 16.71 we have 
unique scalars b^j satisfying 

n—k 

y hi ^ ^ ^ bi^i-\-jhi-\-j^ bi_^i—Yn—k')bi_^i—k 7^ 0. 

j=n—m—k 

m+fc 

By expanding it around P, we get bfi_m-k = Vodi/di^m-k = Vo/fc independent of i. 

Dehne an inhnite matrix A = {aij)i^\n],j& by 

{ m+fc 

fc 

bij+j— - for - m - k < j < -fc, 

Vo 

0 otherwise. 

By using 03.31) for this A, we obtain L{A;x) E (d~^{f). 
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6 .6. Proof of Theorem 15.4L Fix L{x) G and assume that a~^{L) contains 

Q E Ai. By Lemma [7.11 the action of the hnite symmetric group &m via the i?-matrix 
action of ^2.21 preserves the Lax matrix L{x). (See §7]for more discussion of this action.) 
Thus &m ■ Q C a~^{L). By fl2.ip . the &m action induces the permutation action on the 
m quantities ei(g), ... which corresponds exactly to the permutation action on 

Ra- Also, for / G V, these m quantities ei(g), ..., 6^(0') are distinct. So \&m ■ Q\ = 
ml = |a“^(L)|! implying that &m ■ Q = Thus the map a~^{L) -E- {L,Ra) is an 

injection, and the claim follows. 


7. Actions on M 

We study the family of commuting extended affine symmetric group actions on Ai, 
introduced in § 12.21 We also study another family of actions on Ai, which we call the 
snake path actions. The main result in this section is Theorem 17.31 

Throughout we £x / G V. For simplicity, we write Qi := Qi{x) in the rest of 

this section. 

7.1. The behaviour of spectral data under the extended affine symmetric 
group actions. Recall the dehnitions of energy and i?-matrix from 1 12.11 Let A{x) and 
B{x) be two n by n matrices of the shape fl3.ip . with diagonal entries a = (oi,..., a^) 
and b = (6i,..., bn). Dehne the energy by E{A, B) := E{a, b) and the R-matrix to be 
the transformation {A{x), B[x)) i-A {B'[x), A'[x)), where B'{x) and A'{x) have diagonal 
entries b' and a' respectively, and i?(a, b) = (b', a'). The following is proved for example 
in |LP131 Corollary 6.4]. 

Lemma 7.1. Suppose R{A{x), B{x)) = {B'(x), A'{x)). Then we have A{x)B{x) = 
B'{x)A'{x). 

Indeed, Lemma 17.11 and the condition that R is non-trivial uniquely determine the 
R-matrix as a birational transformation. 

Now we reformulate the W xW actions introduced in ^ 12.21 in our current terminology. 
The action of the extended affine symmetric group W is generated by (1 < i < m — 1) 
and IT, where for 1 < i < m — 1, Si acts by 

Si ■ {Qly Q 21 • • • 1 Qm) ' t (Ql, . . . , Qi, ■ ■ ■ ■> Qm) 

where {Qi+i,Q'i) is the R-matrix image of {Qi,Qi+i). The operator tt acts as 

vr : (Ql, Q 2 ,.. . , Qm) ^ {Q 2 , Q 3 , • • • , Qm, P{x)'^QiP{x)-’^). 

The operator e„ acts by moving, via the R-matrix, the last u terms of 

{Qu+l, ...,Qm, P{x)'^QlP{x)-'^, ..., P{x)’^QuP{x)-’^) 


to the hrst u positions, keeping them in order. 
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Recall from 1 12.2l tliat we also have Q = (Qi ,... ,Qn) = {Qij) coordinates on Ai, where 
Qi '■= Qi{x) is an mn' by mn' matrix: 


( ^,1 

0 

0 

X 

\ 

1 

9i,2 

0 

0 


0 



0 


V 0 

0 

1 

Qi^mn' 

/ 


Similarly, the action of W generated by (1 < r < iV — 1) and vr on Q = (Qi, ■ ■ ■ ,Qn) £ 
At is described as follows: for 1 < r < iV — 1, acts by 

Sj. . (Ql, • • • , Qn') ' t (Ql, • • • , Qr+li Qri ■ ■ ■ y Qn) 
where (Q^^Qr+i) R-matrix image of {Qr+i,Qr)- The operator tt acts as 

^ : (Qi, 4,..., Qn) ^ (4, 4, • • •, Qn, P{xr^'QiP{x)-^'^'), 

where P{x) is the mn' by mn' version of P{x) fl3.4p . 

Recall that we denote by &m C W (resp. &n C. W) the hnite symmetric gronp 
generated by Si,..., Sm-i (resp. Si,..., sat-i). Then W is generated by Gm and vr, and 
W is generated by ©at and tt. 

Recall that in § 12.41 we dehne divisors Ou and Au on Cf as follows: 

N u 

Ou = uP- ^ Oi, Au = uP-'^Ai, (7.1) 

i=N—u-\-^ 

for u G Z>o. Let r acts on Sf by (ci,..., cm) eA- (c^f, ci,..., cm-i)- 

Theorem 7.2. (i) The actions of Gm a,nd vr on induce the following trans¬ 

formations on Pic^(C/) X Sf X Rq x Ra-' for r = 1,... ,m — 1, Sr induces 


([P],c,0,Al)^([P],c,0,z/,(Al)), (7.2) 

and 71 induces 

{[D], c, O, A) HA {[D - ^i], r"^(c), O, Vm-ii'm -2 ■ ■ ■ i^i{A)). (7.3) 

(ii) The actions of G^ and tt on 'ip~^{f) induce the following transformations on 
PiA{Cf) X Sf X Rq X Ra: for r = 1,..., N — 1, Sr induces 

{[V],c,0,A)^ {[V], c, h^-riO), A), (7.4) 

and TT induces 

([T>], c, O, A) HA {\p + C>i], r(c), ui--- vn-i{0), A). (7.5) 


Theorems 17.21 will be proved in 1 17.51 The following theorem states that the commnting 
Z™ and Z^ actions on Pic^(C'/) are linearized by the spectral map 0. 

Theorem 7.3. The following diagrams are commutative: 

p-\f) — Pic"(C’/) xSfXRoxRA (7.6) 

eu (—[Au], T““, id, id) 


p-pf) — W{Cf) xSfXRoX Ra, 
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for u = 1,... ,m, and 

— Pic"(C'/) xSfXRoxRA (7.7) 

eu {+[Ou],T'^,id,id) 

— Picnicf) xSfXRoX Ra, 

foru = 1,.. ., iV. Here — o,nd +[0„] respectively act on Pic®(C/) as [D] i-A [D — Ay\ 
and [D] ^ [D + Oy\. 

Proof. We shall show fl7.6p . The proof of fl7.7p is done in the similar way by exchanging 
the rules of the special points for that of {Oj}jg[ 7 v]- 

From fl7.2p and fl7.3p . it follows that a non-trivial part is to show fl2.8p induces the 
action on Pic®(C/) x Ra'- 

{[D],A)^{[V-Au],A). 

Note that for 1 < i < j < m, ■ z/j G &m acts on Ra as 

(Ai,..., Am) eA {Ai , . . . , ^2—15 * * * 5 -^j'l * * * 5 -^m ). (7.8) 

Thus, fl7.3p indicates that vr induces ([P], (Ai,..., Am)) (["D — ^i], (^ 2 ,..., Am, ^i)), 
and vr" induces 


i[P^]i • • • ) P^m)) ([ 7 ^ “ Ay\, {Au+l, . . . , Am, Ai, A 2 , ■ ■ ■ , Am)). 

Further, fl7.2p denotes that does not change a point in Pic®(C/) and acts on Ra as 
a permutation G &m- The rest part (sn • • • Sm-i)(su-i • • • Sm- 2 ) • • • (si • • • Sm-u) of e„ 
changes {A^+i,..., Am, Ai, A 2 ,..., Am) G Ra following (the inverse of) fl7.8p as 

—Li 


• • • 7 ^21 • • • 5 -^m) * ^ • • • 5 -^m-} ^2? • • • 7 -^m) 

^2 ■ ■'^771 —L4 + 1 ——2 / A A A A A \ / A 

HA ••• HA [Ai, . . . , Au-l, Au+1, ■ ■ ■ , Am, Au) HA (^l,-- 


, Ar. 


□ 


We remark that and act on the sets Ra and Rq as the identity transformations. 


7.2. The shift operator. We dehne a shift operator acting on Ad as 

? : (Qi)l<i<m t {P{x)QiP{x) ^)l<j<m- (7-9) 

It is easy to see that it induces an action <;* : C ^ C given by : L{x) ha PL{x)P~^. 

The following result generalizes m Proposition 2.6]. 

Proposition 7.4. The following diagram is commutative: 

fj-^f) Pic<^{Cf) X Sf X Ro X Ra (7.10) 

? {+[P-ON],'r,i>i---C'N-l,id) 

i^-\f) — Pic"(C'/) xSfxRoX Ra, 
where +[P — On] acts on PiO(Of) as [V] [D + P — On]- 

Proof. Suppose Q G '0“^(/) and (j){Q) = {[R],c, O, A). Then (j)og[Q) = {[V'],t{c), z>i • • • z>a7_i(0 ). A), 
where V is some positive divisor of degree g. We set {0[,..., O'j,^) := z/i • • • i>N-i{Oi,... ,On) = 

{On, Oi, ..., On-i)- 
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Recall that g = {gi,g 2 , ■ ■ ■ ,gn)'^ denotes the eigenvector of a{Q). An eigenvector of 
aog{Q) is {g'l, g 2 , ■ ■ ■, g'nV ■= {xgn, gu g 2 , ■ ■ ■, gn-iV■ Dehne h' := g'Jg'^ = gi_i/gn-i = 
hi-i/hn-i] these ratios do not depend on which eigenvector of a o g(^q) we chose. Then 
by dSl]) 

n 

{h[) = (h,_i) - = (A-i -V-{n-i + l)P + Y,0,)- (D„_i -V-P + Or,) 

j=i 

n 

= A-i - Dn-i -{n-i)P + ^ O'. 

j=i+l 

By the nniqneness of V it follows that the divisor V (resp. D[) is eqnal to Dn-i (resp. 
Di_i). Furthermore, V — V P — On, as required. 

□ 


7.3. W X W actions on Lax-matrix. For r = 1,..., — 1, i = 0,..., n' — 1 and 

s = 0,..., m — 1, using the energy in § 12.11 dehne 

Dehne an n by n matrix Bg^r by: 

ys) _ ar+iiq) - ar{q) 


r> _ TT I \ ^ 

■^s,r — ~r y ^ rY^_|_^^Xv7-_|_fc2,r+fc2+l5 


-1(5) 


2=0 ^r-\-ki 

where {Eafi)c,d = ^a,c^b,d for a, b,c,dG Z/nZ, and denotes the identity matrix. 

Lemma 7.5. (i) The action of Sr {r = 1,. .., m — 1) on Ai induces the identity map 

on C, and the action of ti on M induces the adjoint transformation 

IT* : L{x) i-A- Qf^L{x)Qi 

on C. 

(ii) The action of Sr and n on Ai induces adjoint transformations on C, given by 

s* : L(x) !-)■ (Bo,r)~^T(x)Bo,r, 
for r = 0,..., N — 1, and 

TT* : L{x) !-)■ P{x)L{x)P{x)~^. 


Proof, (i) Due to Lemma im we have a o Sr{Q) = o(Q) for 1 < r < m — 1. The induced 
action of vr is 

vr*(gi • ■ ■ QmP{xf) = Q2 ■ ■ • QmP{xfQl, 

and the claim follows. 

(ii) By [LP13( Theorem 6.2], the action of Sr on Ai is given by 

{,Q s) s=l,...,'m ^ ^QsBgj^i (7.11) 

By definition satisfies - £^^+1(^-1)’ £s+m,r — 

The claim follows. 


The second part is easy to see, as the transformation tt just cycles the indices inside 
the Qi fPehnition 12.3p . □ 
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Since similar matrices have the same characteristic polynomial, we obtain the follow¬ 
ing: 

Corollary 7.6. For each q G Ai, ipiq) is invariant under the action of W x W. In 
particular, the coefficients of f{x,y) = fj^q) are conserved quantities for W x W . 


Lemma 7.7. Each c ^ Sj is invariant under actions of Si and Si. The cyclic shift tt 
acts on Sf via 

T ^ : (ci,..., Cm) (c2,..., Cm, Ci), 

and TT acts on Sf via 

T : (ci,..., Cm) {cm, Cl,, Cm-i)- 


Proof. From the definition of the cyclic gronp action fl2.7p . it follows that TT{Q)ij 
and n{Q)ij = qi+ij- Thus, using Lemma 1^751 we see that tt induces 


C£= 

mod M 




1 -^ 


E 

mod M 


Qi,j-l — 




Q—15 


i+j=£ mod M 


and that tt induces q i-A- similarly. 




□ 


7.4. Eigenvector change under conjugation. For a vector v = {vi,V 2 , ■ ■ ■ ,Vn) of 
rational functions on Cf, dehne 

D{v) = (common zeroes of Vi) — (common poles of Vi) + R + On 

where R is the divisor supported on Sp^o = {P, Oi ,..., O^}, which appears in Theorem 
15.11 (See Lemma [6.41 for the explicit formula of R.) The following result is immediate. 


Lemma 7.8. Let f be a rational function on C. Then D{{fvi,...,fvn)) is linearly 
equivalent to D{{vi,... ,Vn)). 


Lemma 7.9. The positive general divisor V = V{L{x)) of degree g associated to L{x) 
is equal to D{{Ai^n, ■ ■., A„,„)). 

Proof. By Theorem 15.11 R = Dn belongs to the common zeroes. Also by Theorem 15.11 
Di,..., Dn have no common points, so they do not contribute to Zi)((Ai ..., A„^„)). □ 

Lemma 7.10. Suppose M = M{x, y) is an xn matrix whose entries are polynomials in 
C[x, y], thought of as rational functions on Cf with poles supported at P. Let D{M ■ v) — 
D{v) = Z1+ — D_, where Dp and D_ are positive divisors. Then restricted to Cf \ P, 
we have that (detM)o — Dp is a positive divisor. Also, is supported at P. 


Proof. Let p G Cf \ P. Then the entries of M are regular at p. We shall show that 
the multiplicity of p in Dp is less than the multiplicity of p as a zero in det M. Since 
d = Mij{.^,y)'>^j, cind Mij is regular at p, it is clear that mn\tpD{Mv) > multpZl(n), 

where multp denotes the multiplicity of a divisor at a point p. We also have 

Since \M\ff[x, y) is regular at p, we have multpil(n) > multp(Mn) —multp(det M). Both 
claims follow. □ 
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Suppose L'{x) = Qi^L{x)Qi and V = 'D{L'{x)). We shall compute "D' — "D up to 
linear equivalence. Note that v = (Ai ..., is an eigenvector of L{xY' — y, and 

L'(x)^ = QjL{x)^, so an eigenvector of L'(x)^ — 2/ is equal to v' = Qjv, where 


Ql = Ql (x) 


/ 

1 

0 

0 

\ 

0 

qi,2 

1 

0 


0 



1 


\ X 

0 

0 

qi,n 

/ 


Extend the vector v into an inhnite vector v by Vi+n = xvi. As in ^ let A = (ajj)ig[„] 
be the inhnite “unfolded” version of L'^{x), satisfying (L^)jj(a:) = (Note 

that A is a matrix of scalars, but h = (... ,n_i,noWi) • • 0^ is a matrix of functions.) 
Then we have 


A - V = yv. (7-12) 

Dehne w = {vi,V 2 , ■ ■ ■ ,Vk+m)^ and w' = (n(, ^ 2 ) • • • claim that there 

exists a (fc + m) X (A; + m) matrix M = M{y) such that w' = M{y) ■ w. Since v' = Qjv, 
we have n' = qiVi + Uj+i (where qi is extended periodically if i > n). We now write 
Vk+m+i in terms of Ui,..., Vk+m using (I7.12p . The matrix A is supported on the m + 1 
diagonals k,k + l,... ,k + m. The matrix A — y is supported on the {k + m + 1) diagonals 
0,1, ■ ■ ■ ,k + m. Thus fl7.12p gives 


ni,fc+iPfc+i + ■ ■ • + ai,fc+m+iPfc+m+i — yxi- 

Also note that ai^k+m+i = 1. So 

/ qi,i 1 0 . 

0 gi,2 10 0 


M{y) = 


0 

0 


0 

0 


0 


0 


qi,k+i 


(7.13) 


\ 


\ y 0 ■ ■ ■ -ai,fc+i • • ■ 

Lemma 7.11. With the above conventions, we have 

V r^V + Al-P. 


Qk+m —1 1 

Uljfc+m—1 qi,k+m, 0,1,k+m ) 


Proof. We shall show that D{y') — D{v) = Ai — P. This suffices by Lemmas 17.81 and 17.91 
From Lemma 17.101 we have D[v') — D[v) = — D_ where D± are positive divisors, 

with Z1+ when restricted to Cf\P supported on {(detQ^)(a:) = 0}, and supported 
at P. Since D{v') is a positive divisor of degree g by the construction of L'(x)^, we have 
that Z1+ and have the same degree. 

We now calculate D_. By Theorem 15.11 we have that multp(nj) = C — i for some 
integer C (we use the convention that negative multiplicity is a pole), and this formula 
still holds for the inhnite vector (... ,V-i,vo,vi ,...). Since n' = qiVi + Ui+i, we have 
multp(n') = C — {i + 1). Thus D_ = P. It follows that Z1+ is a single point in Cf\P. 

We shall show that Z1+ must be supported on Ai. By Theorem 15.11 the common 
zeros of Vi and Vj for any i ^ j are only Dn except for the points in Sp^o- Thus we 
have D{w) = Dn + R' where R' is supported on Sp^o- But {(detQ^)(x) = 0} does not 
















GENERALIZED DISCRETE TODA LATTICES 


35 


intersect Sp^o- Using Lemma IT. 101 now applied to w' = M{y)w, we see that -D+ must be 
supported on {det M{y) = 0}. 

By Lemma [7.121 below, we conclude that Dp is a multiple of Ai. Since D{v') — D{v) 
is a divisor of degree 0, we must have 

D{v') - D{v) = Ai-P, 

as required. □ 

Lemma 7.12. The intersection of {detQ^){x) = 0 and {det M{y)) = 0 is the single 
point Ai. 


Proof. We check that det M{y) = ±y. It is clear that 

/ Qi,k+i 1 

Ql,k+2 1 

det M{y) = (-l)”"+^|/+gi,i • • ■ qi^k det 




Ql,k+m—l 1 

V ®l,fc+2 ' ' ' Ojl,k+m—l Ql,k+m 0‘l,k+m/ 

Write Ui {i = 1,..., m) for the m columns of the above m by m matrix. We show that 
they are linearly dependent as 

Q'l,fc+l (^2 Q'l,fc+2(^3 ■ ■ ■ igJ'm—l gi,k+m—l'^rri) ' ' ' )) ^5 

which reduces to 


~ + ?fc+l®l,fc+2 ~ Q'l,fc+l5l,fc+20l,A:+3 + • • • ('^•14) 

“1“ ( 1) <?fe+l<?fe+2 ■ ■ ■ gk+m—lOil,k+m A ( 1) Q'fc+lQ'fc+2 ‘ ‘ ' Qk+m ~ 0. 

By the dehnition of A, we have 

J/.i j = k + l,...,n, 

aij = < 

j =n + l,...,n + k, 

where L = {lij)i^[n],jez is the inhnite unfolded version of L{x). On the network G', lk+i,i 
is the weight generating function of paths from the source k + i to the sink k + 1. Note 
that each weight in lk+i,i is of degree m + 1 — i in q-s. We divide lk+i,i into two parts: 

lk+i,l = lk+i,l A 

where I'k+i^i is the sum of the weights with qi,k+i, and I'k+i^i is the sum of the weights 
without qi^k+i- We claim that there is one-to-one correspondence between the paths con¬ 
tributing to I'kpi^i and the paths contributing to Precisely, we have qi^k+d'k+i+i,! = 

lk+i,i- Combining with = 0 and 4+^1 = qi,k+m, we obtain fl7.14p . □ 

7 . 5 . Proof of Theorem [731 (i) Due to the dehnition of Sr, Lemma I775f il and Lemma 1771 
for Q = (Qi, ..., Qm) e with 0(Q) = ([D], c, O, (Ai, ..., A^)), we have 

0 O Sr{Ql, . . . , Qm) = ([7*], c, O, (Ai, . . . , Ar-I, Ar{Q'), Ay-pi^Q'), ^^+2, • • • , ^m)); 

where Q' = Sr{Q). From fl2.ip we see that er{q') = er+i{q) and er+i{q') = 6 ^( 9 ), and fl7.2p 
follows. 

It is easy to see that for L{x) = a{Q) we have y o ■n[L{x)) = ([D'], r“^(c), O'), where 
D' is some divisor. By Lemma [7. Ill we have [D] = \D' + Ai — P]. The claim follows. 
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(ii) First we prove fl7.4p . Let g = {gi,... ,gnY be the eigenvector of L{x) G V' and 

set , g'^^ := Define h[ := g'i/g'n for i = 1,..., n — 1. Then we have 

/ _ J gr+ki - j = r + ki mod n 


g^ = 


gj 


otherwise. 


We mnst show that (h' )oo = (hj)oo for r = 1,..., A^ — 1 and j = 1,..., n — 1. It is enongh 
to consider the case of r = 1 and j = 1. Dne to Theorem 15.11 and fl5.ll) . there are positive 
divisors D, D'of degree snch that (hi)oo = {n — l)P + V, (h'^)oo = {n — l)P + V'. Since 
h[ = hi- Kf'*h 2 and (hi)^ > (K®h 2 )oo = (^ 2 ) 00 , we get (h'Joo < (hi)oo- Thns it follows 
that V' = V. From fl7.1ip and Lemma 17.71 it follows that both A G Ra and c G are 
not changed by Sr- From the facts that ajuj^i^rig) is the prodnct of all diagonal elements 
of Qr, and that the i?-matrix action changes {Qr,Qr+i) to {Qr+i,Q'r), h follows that Sr 
exchanges the {N — r)-th and the (iV —r + l)-th elements of O G Ro- Then we get fl7.4p . 

Since tt indnces the shift of L{x), 07.51) follows from Proposition 17.41 


7.6. Snake path actions. Recall that M = gcd(?7,, k + m). The snake path actions are 
torus actions on At, considered in |LP13j . Recall from li l4.2l that a snake path is a closed 
path on G that turns at every vertex. Thus it alternates between going up or going 
right. For 1 < s < M and t G C*, the action Tg := Ts{t) on At is given by 


tqij 

if j = s — i + 1 mod M 



ii j = s — i mod M 

(7,16) 

Qij 

otherwise. 



Informally, Ts{t) multiplies all left turns on a snake path by t, and all right turns by 1/t. 


Lemma 7.13. The induced action T* (1 < s < M) on C is given by 

T; : L{x) ^ Dg{t)L{x)Dg{t)-\ 

where 

T/ 7 \ 1 1^ if i = s mod M 

Ds{t) = diag(ddi<i<n, di = < 

I 1 otherwise . 

Proof. For simplicity we write Dg := Dg{t) and P := P{x). We rewrite DgL{x)D~^ as 
DgL{x)D:^ = i^{p-^+^DgP^-^)Q,{p-^DgPY^ ■ {P-^DgP^)P'^D:\ (7.16) 


This is equal to a{Tg{Q)) since (P ^DgP^) ■ P^D^ ^ = P^. 


□ 


Proposition 7.14. We have the following commutative diagram: 


At At 



At -- At. 

<r-“oTsO<r“ 
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Proof. It is enough to prove ((<7*) “ o T* o (<7*)“) o e* = e* o T* acting on £. For L(x) = 
Qi ■ ■ ■ QmP^ £ P) we have 

L{x)y^Qu+i---QmP''Qi---Qu 

... Q^pkQ^... 

m 

= n {P~"^^DsP"~^)Qi{P~"DsPY^ ■ {P~'^D,P^)P'^Df^ 

i=u-\-l 

u 

■ J](p-*+ip,P*-')Q,(p-*P,PO-' 

i=l 

where Q[ = {P~''^^DsP''~^)Qi{P~'‘‘DsP'^)~^. On the other hand we have 

L(x) A D.L(x)D:' QIP>= 

Aq'„+i...q'„p‘q',---Q7 

Thus the claim is obtained. □ 


From Lemma [6.51 and [Y.13l we obtain 


Proposition 7.15. The following diagram is commutative: 

V'-H/) — Pic"(0/) xSfxRoxRA 

Ts {id, ts, id, id) 

fj-Hf) Pic"(C'/) xSfxRoX Ra, 

where tg ■= ts(t) acts on Sf by 

{ t Cl if s = i 

j Cl if s = i + 1 
Cl otherwise. 

Thus the snake path actions Ts{t) act transitively on iSj. 


(7.17) 


(7.18) 


8. Theta function solution to initial value problem 


8.1. Riemann theta function. Fix / G V. We £x a universal cover of Cf and Pq G Oj, 
and dehne the Abel-Jacobi map l by 








where ooi,... ^Ug is a basis of holomorphic differentials on Cf. We also write l for the 
induced map Div°(0/) —?■ C^. Let O be the period matrix oi Cf, and write 0(z) : = 
0(z; O) for the Riemann theta function: 


0(z; O) = exp (tt-v/— Im ■ (Om + 2z)) , z G O. 











38 


REI INOUE, THOMAS LAM, AND PAVLO PYLYAVSKYY 


It is known to satisfy the quasi-periodicity: 

0(z -I- m -I- nf2) = exp (—7r\/—1 n • (fin -|- 2z)^ G)(z), 

for m, n G Z^. 

Let T>o G Div^“^(C'/) be the Riemann characteristic, that is, 21^0 is linearly equivalent 
to the canonical divisor Then the well known properties of the Riemann theta 

function gives 


Lemma 8.1. For any point Y G Cf, the function 

X ^ 0(i(x -y)) 

is a section of the line bundle 0{Vo + Y) which has zeroes exactly at Vq and at Y, and 
no poles. For any positive divisor V of degree g, the function 

X ^ e{L{x -V + Vo)) 

is a section of the line bundle OifD) which has zeroes exactly at V and no poles. 

8.2. The inverse of 0. For a positive divisor T) of degree g, dehne the function fji {i G Z) 
on Cf by 

, ^ 0(/3 - 6(a-.) + i{x - o^)) nL+i mx - o,)) 

’ 0(/3 + t(X-O„))0(6(X-P))-- 

where fd = l{Vq + On — V). 

Lemma 8.2. fji{X) is a meromorphic function on Cf. When P is the positive divisor 
defined by Definition \5.S\. we have 

n 

m = {K) = V,+ 0,-V-{n- i)P. (8.1) 

j=i+l 

Thus 'ifi{X) is egual to hi{X) up to a scalar. 


Proof. To check that '4)i{X) is a meromorphic function on Cf we use the functional 
equation for 0, and note that the multiset of points (with signs and multiplicities) that 
appear in the numerator is equal to the same for the denominator. This multiset in 
additive notation \s fd — {n — i)P + {n — i-\- 1)X — On- For the second statement, we note 
that 'ifi{X) has zeroes at Oj+i, Oj+ 2 ,... ,On, a pole of order n — i a.t P, and also poles 
at V. By 05.111 . we have = (hi). □ 


Now we study the inverse of cf. We focus on the Z™' action. For t = (ti,t 2 ,..., tm) G 
Z™, let g* := (gjjji G 'ip~^{f) be a conhguration at time t. For P{x) := a{q^), let 
g* = {gl,... be its eigenvector and set h\ := gl/gh- Let 4){q^) be 

(P , (c;^, • • • , Cj^), (Ai,..., Am), (Oi,..., On)), 


where = P — some positive divisor P of degree g. Dehne e^- : = 

(1,..., 1, 0,..., 0) for j = 0,..., m — 1, and recursively set = (1,..., 1) -|- . Dehne 


J m-3 

functions 6\ and T* on C/ by 


e\{X) = e{(d^ + L{X-On-On-f)), 
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where = t(Po + On — V^), and by 




t+Sj 




«‘+'“-(.y)CSA)' 

By Lemma [8.21 there is a constant depending on C/, such that we have 




Lemma 8.3. (Cf. Lemma 3.1 mi) We have 






D,* 




t+e,' aj-Af-l 

“i-1 “i-W 


( 8 . 2 ) 

(8.3) 


where d\ is a eoefficient of the leading term of h\{X) around X = P. 


Proof. The first equality of fl8.2p follows from 6\{Oi) = 6\_^{P) and O* = Oi_Ar. By 
Lemma [7.51 we have = Q* which gives 

= B+e, (C^-l(X) + (X))hr- (X) 

Then, from fl8.ip the second equalities of 08.211 and 08.3p follow. 


Lemma 8.4. VTe have 


d, 


,t+ew 


d 


i+l 


i + j = i mod M, 


□ 


(8.4) 


where we extend cj to i E Z by setting c\ ■= c\ if i = £ mod M for 1 < £ < M. 


Proof. It is enough to show the case of j = 0, 1 < i < M. Then the other cases follow 
from Lemma [7.71 

We will show that the coefficient of the leading term in at X = P is equal to Cj. 
To obtain the ratio we may choose to compute maximal minors of L{x) — y with 
respect to any row. Let us choose row i + l. Then by Theorem 16.21 we are counting 
families of paths that do not start at source i + l and do not end at sink i (for gi) or 
sink i + l (for Qi+i). 

Since we are computing at X = P, the only terms that contribute are the ones lying on 
the upper hull of the Newton polygon, that is, the edge with slope —n/{m + k). In terms 
of path families this means that we only take paths that are as close to a snake path 
as possible. In the case of this means that we take a subset of closed snake paths. 
Each of them picks up no weight, and thus overall we just get a constant. In the case of 
gt one of the closed staircase paths skips one term, and thus we get a long path starting 
at source i, winding around the torus and ending at sink i + l. Such path essentially 
by dehnition picks up weight c*. The other snake paths may or may not appear, thus 
creating only a constant factor in front. □ 


Note that fl8.4p is compatible with the snake path actions in Lemma 17.131 
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Theorem 8.5. When N = 1, the inverse map of cj) is given by 


Wi = fc 


■ Qj ■ 


where 


aj = x{Aj) ” ■ exp 


o: 


iP)dZ?iP) 


n 


h-i(P- A,) 


(8.5) 


and b ^ is determined by a,b G TA such that a + bQ := i{On)- 
Proof. Let p* ^ = q]^iq]^i_i ■ ■ ■ Q'j j_jv+i- By Lemma [H3] we get the equations 


Pii 


j ' •,Tr^+®7 / T~>\ »Tr^+®7 / T~>\ J^+®7 


^rip)G ^r-?iP)G-i 


We show that aj does not depend on t. Consider the product 


^(^j ) Pj,iPj,i+N ’ ' ' Pj,i+(n'—l)Ni 


which is equal to 


a 


t\n 




Uj, 


nt+f^j-l ( TD\ ( TD\ 

i+n-N '^i+n-N W ) ) 

0 / 


nt+Gj 


Since On is equivalent to 0 in Pic‘^(C'/), there exists a, fe G such that a + hVL = i{On). 
Due to the quasi-periodicity of Q{z), we have 

9l^{P) = exp [-27iV^b • (/3‘ - L{On-^ - Oi) + bn/2)] • 0j(P). (8.6) 

By using fl8.4p and fl8.6p we obtain aj independent of t as 


_ N_ 

aj = fc "af. 


Hence we get 


Pli = f^ 


N 

■"af 


N 

n 


»P(P)»ZT'{P) 


t+Gj-l . 


1=1 


When TV = 1, it is nothing but the claim. 

8.3. Conditional solution for iV > 1. When > 1, by factorizing fl8.7p we obtain 


(8.7) 

□ 


Qj,^ = fc^- Ipi • a,- • C+i-i^+i- 


0 j+"^(P) e^iti~\P) 


Here satishes 


N 


^=1 


i+l 


e7"‘-\P)e‘*‘’(P) 


= 1 , 


(8.8) 


(8.9) 


which denotes that = 7j,i+Af- 

Lemma 8.6. For i such that i = r mod N we have 


i=i 


n 

i=i 


-1 




P, 


exp ( 27r\/—1(6' — b —) ■ l(P — OP 
\ n' 


( 8 . 10 ) 


Here b' is given by a', b' G Z^ such that a' + b'Q := L{Am + Ok). 
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Proof. From fl3.6p and 


it follows that 

m n'—l 


^ j nn 

^jr' = l / j = l p=0 

for r = 1,..., iV and i = r mod N. 


, ^^' oh,N-iiP)on]NiP) 

^j+i+pN—l J_J_ nt+Bm {]D\ Qt I pi ’ 
j=0 ) ^i+jN\-^ ) 


( 8 . 11 ) 


First we show that the second factor of 08.1111 is equal to . It is enough to prove 
in the case of f = 1, where the second factor is 11 ^= 0 ^ (^i+pN • • • Cm+pN)- When n’ = k', 
M is a divisor of m and we have ci+pN ■ ■ -Cm+pN = When n' ^ k', define an 

automorphism u of A/"' := { 0 ,..., n' — 1 } hy u : p ^ p+n’ — k’ mod n'. Since n' and k' are 
coprime, a set {z^^(l) | i. G A/"'} coincides with M'. Thus, from M = gcd(n, k+m) and m+ 
pN = i>{p)N mod M, it follows that the product Y\^=o i.^i+pN''' Cm+pN) is reordered 
to be Ylj^=i ^ 3 - Further, since M is a divisor of mn', we obtain the claim. Next, we can 
show that the third factor of 08.111) is equal to exp (27r\/—1(—fe'n' + hk') ■ l{P — Or)) , 
by using 08.6p and 


0i+e.(p) ^ [_2Tr^b' • (/3* - L{On-^+k - Oi) + &'fl/2)] ■ 


Here b' is given by a', b' e such that a' + b'kl := L{Am + Ok). Finally from 08.lip we 
obtain 



and the claim follows. 


= (Tr ■ ttj ” • exp —l{n'b' — bk') ■ l{P 
i=i 


Or)), 


□ 


Unfortunately, we have not been able to obtain from 08.9p and O8.10p . Nevertheless, 
if we assume that is constant on Pic^(C/), then we obtain the following conditional 
result. 


Proposition 8.7. Suppose that is a constant function on Pic^(Uj). Then we have 


lj,i = o-p "‘ ■ Oj • exp 
i=i 

In particular, the inverse of 0 is given by 


(— y {n'P — k'h) ■ l{P — Or) 

\ n’m 


( 8 . 12 ) 


^ 0r^(P)Ci^-(P) 

q^. = Q.ayOi- ^t+e. 


e^^-^p) cr(^) 


where aj is defined by 08.5p . and 

m 

Q = fc^ 'YlxiAj 


1 

Tim 


= (- 1 ) 


i=i 

\n'+iyiOr)"' 


x(Or 


\k’ 


1 

n'm 


• exp 




nm 


{{nb'-kb)-L{P-Oi)-b-L{AJ) 


i = r mod N. 


Proof. To have 08.8p compatible with the snake path action 07.151) and 07.181) . has to 
be constant on Sf. So is regarded as a function on Ra x Rq. If 7j,i is not constant 
on Ra, 08.1 op implies that 7 j,j depends on af^, but this contradicts 08.Op . Thus 7 ^,* is 
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constant on and we have 08.121) which fnlhlls 08.9p . By snbstitnting 08.12p in 08■8p 
and nsing 03.lip , we obtain the hnal claim. □ 


9. Octahedron recurrence 

We shall prove that the fnnction 9\{P) satishes a family of octahedron recnrrences 
Spe| , as a specialization of Fay’s trisecant identity. We follow the dehnitions in ^ In 
the following we assnme N = gcd(?7,, k) = 1 and write O for the nniqne special point Oi 
over (0, 0) E Cf. 

9.1. Fay’s trisecant identity. We introdnce Fay’s trisecant identity in onr setting. For 
a, /5 G M®, we dehne a generalization of the Riemann theta fnnction: 

©[a,/S](z) := exp (vr-v/— 1/3 • (/3f2 + 2z + 2a) ■ ©(z + f2/3 + a). (9.1) 

We call [a,/3] a half period when a, (3 E (Z/2)^. Fnrthermore, a half period [a,/3] is odd 
when Q[a,(3]{z) is an odd fnnction of z, that is, Q[a, (3]{—z) = —Q[a,(3]{z). We note 
that the Riemann theta fnnction itself is an even fnnction: ©(z) = ©(—z). It is easy to 
check that a half period [a, /?] is odd if and only if 4a • /? = 1 mod 2. 

Theorem 9.1. |Fa//| , §/// For four points Pi, P 2 , Pz, on the universal cover ofCf, 
z E C®, and an odd half period [a,/ 3 ], the formula 

Q{z + i{Pz - Pi)) Q{z + i{P, - P 2 )) Q[a,P]{i{P 2 - P 3 )) Q[a,(3]{i{P, - Pi)) 

+ ©(^ + i{Pz - P 2 )) Q{z + 6(P4 - Pi)) ©[«,/9](i(Pi - P 3 )) 0[a,/5](^(P2 - P 4 )) 

= ©(z + i{Pz + P4 - Pi - P2)) ©(z) ©[a,/3](.(P4 - P3)) ©[a,/3](.(P2 - Pi)). 

holds. 


9.2. m = 2 case. When m = 2, the vertical actions [u = 1, 2) are written as difference 
eqnations expressed as 

92,i^l,i-k — ^ ( 9 - 2 ) 

'?2,i+l + ~ '?l)i+l + 92^1 ^ (9-3) 

<llr-k = <ltr (J = l,2). (9.4) 

For simplicity we write 9\ for 9\{P) = ©(/3* + i((n — i — 1){0 — P))). By constrnction 
the theta fnnction solntion of ^ (Theorem [53]), 




(9.5) 


satishes fl9.2l) - fl9.4p . 

Theorem 9.2. For any t E 7? and i E "L, the dj satisfy an octahedron recurrence 
relation, 


0t+e2 0t+2ei _ 0t+2ei 0t+e2 _ ^ ^t+ei+e2^i+ei 


i+1 


i+1 

Here c is a constant given by 

_ 0(po + ^(^1 ~ ^2))©(Po + I'iO — P)) 
0(Po + /'{O — A2))Q{po + i{Ai — P)) ’ 
where po is a zero of the Riemann theta function: Q(po) = 0. 


(9.6) 
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Proof. By setting {Pi, P 2 , P 3 , P 4 ) = {A 2 , 0 , P, Ai) and using fl9.ip . we obtain: 

T. e(. + i{p - A,)) e{z + -o)) + n e{z + z(p-o))e(z + 

= T3 e{z + i{p PA 1 -A 2 - O)) Q{z) ^ ^ 

where p := PL(3 + a, and 

Ti := Q{p + l{0 — P)) Q{p + i{Ai — A2)), 

T 2 := e4-^/^/3•dA2-Al) 0 ^^ ^ 0 ^^ ^ ^g_g^ 

T3 := Q{p + i{Ai - P)) Q{p + 6(0 - A2)). 

By setting z = (]^ + L{{n — i — 1)(0 — P) + 2Ai), fl9.7p turns out to be 

Ti 0 ^+"^ + T 2 elix = T3 ^. ( 9 . 9 ) 


Lemma 9.3. We have 

T 2 _ Qi 
T3 02 

Proof. For G iJ~^{f), take to £ such that G is a zero of the Riemann theta 
function, which is always possible by choosing (f appropriately. We write —po for such 
. Then we have 9{f_i = 0(— Po) = 0, and q{°^ = q 2 °~^^ = 0. From fl9.2p and fl9.3p . we 
obtain 


„to-ei Ao-ei _ to to 


io-ei _ to 
Hl,n-k-l — H2,n-1- 


(9.10) 


On the other hand, when z = Pq + i{A 2 — P), fl9.7p becomes 

T 2 0(po + ^(^2 ~ O)) 0(po + 6 (Ri — P)) = T 3 0(po + k{Ai — O)) 0(po + tz{A 2 — P)). 

It is rewritten as (using that Q{z) is even function) 

rri flio+ei nto+S2—ei „to—ei 

-‘-2 _ On-2 ^n-1 _ H2,n-1 

F “ eWWF “ F' £7’ 

where we use fl9.5p to get the last equality. It follows from fl9.10p that q^fn-i!Q au-i — 1; 
and we obtain the claim. □ 


We continue the proof of the theorem. By setting z = pq + l{0 — P) in fl9.7p . we obtain 
0(Po + t'{0 — ^ 2 )) 0(Po + ^(^1 ~ -P)) = ^3 0(Po + ^(^1 ~ ^ 2 )) 0(Po + i'{0 — P)). 
Using this and the above lemma, fl9.9p is shown to be 

c()f+« 95, + a, 9 ?=i or = a, OJ, Of', □ 

Conversely, we have the following. 

Proposition 9.4. Suppose that 9\ satisfy 09.61) . Then gC defined by 09.51) satisfies 

(ra-dOD. 

Proof. It is very easy to see that 09.51) satisfies 09.21) and 09.4p . We check 09.3p . We 
consider a ratio 

y2,i+i yi,i+i 
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By substituting 09.51) in /* we obtain 
/• = 


a c - a 

a2Ci+2^^ OlCi+i^eT^ 


, _ . 9T^^^9J]Aa2.9T9j^\-a,9j^^9T.y 

a2Cj+l 2ei ei+ea -fllCj ej .ei+ej ^+l V ^ ^ i ^ 1/ 

^i—1 '~'i '^i — 1 


^ei ^i-1 


(9.11) 


where we omit the superscripts t for simplicity. At the second equality we have canceled 
all the Ci using c®“ = €*+„. Due to 09.61) . we obtain f-=l. Thus, using 09.4p . we obtain 

ra . □ 

9.3. General m case. The vertical actions (n = 1,..., m) are expressed as matrix 
equations: 

g.+e. . , , g^e. ^ g<_^^ , , . . . . QIP(:,)-\ (9,12) 

Among the family of difference equations we will use the following ones later: 


n-'T"- n 4. If? 


i 


(9.13) 


4=1 4=1 

for n = 1, • • • , m. Here we dehne 

(j + n,i) j = 

{j + u — m^i — k) j = m — M + 1 ,..., m. 

Theorem 9.5. For any t G i G Z and 1 < p < r < m, the 9 \ satisfy an octahedron 
recurrence relation, 


s{j,i) = 


A 

^p,r 


'9T 


+ ttp 9, 


‘P '^i+l 


'97 


_ rtl- 


^97 


(9,15) 


Here 6p,r is a constant given by 


(9.16) 


_ 0 (po + i(Ap — Ar)) 0 (po + <-(0 — P)) 

0 (po + yO — Ar))Q{po + i{Ap — P)) 
and po is a zero of the Riemann theta function: 0(po) = 0. 

Proof. The proof is similar to the m = 2 case. We explain the outline. In the case 
(Pi, P2, P3, P4) = {Ap, O, P, Ap) of Theorem 19.11 we obtain 
Ti 0(z + t(P — Ap)) 0(z + t(Ap — O)) + T2 Q{z + t(P — O)) 0(z + i(Ap — Ap)) 

= T 3 0 ( 2 : + i(P + Ap — Ap — O)) Q{z) 

where Ti, T 2 and T 3 are given by 09.8p . but we replace Ai (resp. A 2 ) with Ap (resp. Ap). 
By setting z = A + c((n — i — 1)(0 — F) + Ap + Ap-i) at 09.16p . we get 

rji /qi+ep —l+Gr — 1 /qi+Gp+Gr . rji /qt+ep+Gr — 1/qi+Gp —1 H-Gr rji /q^+Gp —1 H-Gr ^t+Gp+Gr — 1 

^l^z+1 -Js^z+l 

We take to G Zi™' and dehne po := —A° G in the same manner as in the proof of 
Lemma [231 Then = 0 holds for u = l,...,m. From 09.131) (resp. 09.141) 1 of 

i = n — 1 (resp. i = n) and n = p — 1 or r — 1 , it follows that 
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On the other hand, by setting z = po + i{Ar — P) at fl9.16p . we obtain 

rji QtQ-\-Gr—G ^—1 1 

-^2 ^n—2 ^n—1 _ ^r,n—l 

Til ntQ-\-er—er—i n^o-\-Gp—ep—i to—ep—i ^ 

^3 ^n-2 ^n-1 %,n-l 

From the above two relations, T 2 /T 3 = Op/a^ follows. Finally, by setting z = po + l{0 — P) 
at fl9.16p we obtain the formula of 6p^r- D 

The following conjecture extends Proposition 19.41 

Conjecture 9.6. Suppose 6 *- satisfy 09.151) . Then gO dehned via 09.51) satisfy all the 
difference equations 09.12p . 

We have checked the conjecture for m < 3. 


10. The transposed network 

10.1. Transposed Lcix matrix and spectral curve. Corresponding to Q of 02.3p . we 

also identify q E M. with an iV-tuple of n'm x n'm matrices Q := (Qi(a:))jg[Ar] as 1 17.11 
The matrices Qi{x) are given in terms of the by 

QN+i-i{x) := P{x) + diag[gm,j,... ) Qm,i+ki ■ ■ ■ 1 qi,i+ki ■ ■ ■ i qm,i—k) • • • ; qi,i—k\) 

where P{x) is the mn' x ran' matrix: 

0 X \ 

0 0 

-.0 
1 0 / 

We dehne L{x) := Qi{x) ■ ■ ■ Qn{x)P{x)'‘""^, which is another Lax matrix. 

Also dehne a map p : Ai ^ C[x, y] given as a composition, 

Q I—)■ L{x) I—)■ det(L(x) — y). 

Consequently, for g G Ad we have two affine plane curves C^{q) and C^^^^in C^, given by 
the zeros of p{q) and piq) respectively. The proof of Proposition 110.ll is similar to that 
of Proposition 13.31 


P{x) : = 


/o 
1 

0 
V 0 


0 

0 






Proposition 10.1. The Newton polygon N^ip^q)) is the triangle with vertices {0,mn'), 
{mk',0) o,nd {mk' + N,0), where the lower hull (resp. upper hull) consists of one edge 
with vertices {mk',0) ( 0 ,mn') (resp. (mk' + NA) to {0,mn')). 


Lemma 10.2. The affine transformation 


HA 


k'{k + m] 
—n'k 


—k' (i — k'k')/n' 
n' k' 


sends integer points of Nipf^q)) into integer points of Nipp^q)). 


( 10 . 1 ) 


Proof. It is easy to see that it sends the vertices correctly. By the dehnition of k' we 
know (1 — k'k')/n' is an integer. Thus this transformation sends integer points to integer 
points. So does the inverse, as the determinant of the matrix involved is —1. □ 





















46 


REI INOUE, THOMAS LAM, AND PAVLO PYLYAVSKYY 


Example 10.3. Let {n,m,k) = (6,3,4). The two Newton polygons N{'ijj{q)) and 
N{'ip{q)) are illustrated in Figure [TTl Here i labels the horizontal axis and j labels 
the vertical axis. The dots of the same color show integer points inside the Newton 




Figure 11. 


polygon that get sent to each other by the transformation fllU.ll) . The formula for the 
transformation in this case is 



Proposition 10.4. For q G Ai, the polynomials 'ip{q) and 'ip{q) coincide up to the 
monomial transformation induced by fllO.ip . The signs of the new monomials are derived 
from the rule given in Theorem AS.^ the sign of x°‘y^ is 

See § 110.31 for the proof. 

10.2. Special points on the transposed curve. We write f{x, y) for the polynomial 
obtained from the hxed polynomial f{x,y) in ^3.31 via the affine transformation fllO.ip . 
Let Cjhe the smooth compactihcation of the affine plane curve Cj given by f{x, y) = 0. 
As for Cf (Lemma 13.8 p . Cj has a unique point P lying over oo. Due to this fact and 
Proposition 110.41 the two curves Cf and Cj are isomorphic. Let t : Cf ^ Cf he the 
birational isomorphism, which is given by 

(a;,i/)^(l/"V',x^/) 

when {x,y) G C^n (C*)^. Besides P, on Cj we have special points Or = ((—l)”^”'Vr, 0) 
for r G [N], and Ai for i G [m], where near Ai there is a local coordinate u such that 

We see that r(P) = P, r(Aj) = A^ and T{Or) = Or- 

10.3. Proof of Proposition [T074l The terms contributing to a particular coefficient of 
the spectral curve are described by Theorem 13.2[ We shall exhibit a bijection between 
the terms of the coefficient associated with the lattice point (i,j) inside the Newton 
polygon N{fip{q)) and the terms contributing to the coefficient associated with the image 
of that point in N{'ip{q)) under the affine transformation fllO.ll) . 
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An underway path in the network G is the mirror symmetric version of a highway path. 
Thus the weights of an underway path are given by Figure [5] with 0 and qij swapped. 
The crucial observation is that families of highway paths on the network associated 
with Q are families of underway paths on the original network, parsed in the opposite 
direction. This is because by the construction of the transpose map between Q and Q, 
the corresponding toric networks are the same but are viewed from opposite sides of the 
torus (inside vs outside). 

Now, assume we have a closed family of highway paths contributing to one of the 
coefficients of ipiq). Simply complement all edges that ended up on our family of closed 
highway paths inside the set of all edges of the network. We claim that the result can 
be parsed as the desired family of closed underway paths. 

Indeed, the original family can be viewed as a number of horizontal “steps through” 
picking up a weight at some node, connected by intervals of staircase paths that do 
not pick up any weight. We can interpret our procedure as complementing used intervals 
of staircase paths, that is, making them not used, and vice versa. As a result, locally 
around each weight that was picked up the new path will look like what is shown in 
Figure dH Thus, it will be an underway path, and it will pick up exactly such qij. In 
other words, the weight of the original highway family is the same as the weight of this 
new underway family. 


Qij 


qij 


Figure 12. 

It is also easy to see that the new underway family is closed. This completes the proof. 

Example 10.5. Let (a, b, c, h) = (3, 2, 3, 2) as in Example 12.21 Figure [13] shows an 
example of a family of paths contributing to the purple term of the spectral curve as 
marked in Figure [TTJ The hrst step takes the complement of edges of this family inside 
the set of all edges of this network. The second step does not change the network or the 
paths, it just changes the point of view and reverses paths’ directions. 

11. Relation to the dimer model 

In this section, we give the explicit relation between the i?-matrix dynamics on our 
toric network and cluster transformations on the honeycomb dimer on a torus. See EH 
for background on the dimer model. 

11.1. Cluster transformations on the honeycomb dimer. The calculation in this 
section is the dimer analogue of the highway network computation of the geometric 
R-matrix (cf. [LP131 Theorem 6.2]), which was explained earlier in 1 14.51 (Figure [9]). 
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Fix a positive integer L, and consider a honeycomb bipartite graph on a cylinder as 
in Figure [m where {i G Z/LZ) are the weights of the faces in three cyclically 

consecutive rows. We write a := (cii, A, 7j)igz/Lz- 



Figure 14. Honeycomb dimer on a cylinder 

For the weights a of the honeycomb, we define a transformation of a in a following 
way: first we split the ai-face into two, by inserting a digon of weight —1 as the top of 
Figure 1151 We thank R. Kenyon for explaining this operation to us. 

Set the weights of the new two faces to be —c and where c is a nonzero parameter 
which will be determined. Let D be the quiver dual to the bipartite graph, drawn in blue 
in the figure. The weights of faces are to be regarded as coefficient variables associated 
to each vertex of D. With 1,... ,L,a and b, we assign the vertices of D in the middle 
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row, as depicted. Next, we apply the cluster mutations /ii, /i 2 , ■ ■ ■, to {D, a) in order, 
recursively defining oji (i = 1 ,..., L) by 


Oi 


oji —, oJi ai(l + oJi-i). 


The condition that the digon’s weight is again —1 after the L mutations gives an equation 
for c as —c(l + ujl) = —1. By solving it, c is determined to be 

1 - nf=i 


c = 


and oJi := 


1 -r-rt ’ 

, (Tl) is obtained as 


..., ai) 


X]t=0 ns=l 


1 - nti 


a.. 


( 11 . 1 ) 


Definition 11.1. Let Ra be the transformation of a given by 

(oj,/3i, 7*) i-A-(ci;j_\(l + (Ui), /3i(l + 7j(l + a;j^) ^) , (11-2) 

which is induced by the sequence of mutations ■ ■ ■ /^i, 

Due to the expression of oJi, we see that Ra does not depend on which aj-face we split 
at the first stage. 

Remark 11.2. Note that our derivation of Ra is not entirely a cluster algebra compu¬ 
tation since we began with the “digon insertion” operation, which does not have a clear 
cluster algebra interpretation. In an upcoming work [ILPj . we plan to further clarify the 
cluster nature of the geometric i?-matrix. 

11.2. Relation with the (n, m, /i;)-network. We consider the honeycomb bipartite 
graph on a torus as in Figure [HI We assign each face (resp. edge) with a weight Xji 
(resp. qji or 1), satisfying the periodicity conditions Xj^i+n = and Xm+j,i = 

(resp. qj^i+n = <ij,i and qm+j,i = In the figure we omit weights that are equal to 

1. The Xji and qji are related by 

{j = Xm,i = , (11.3) 

hence the product of all the Xji is equal to 1. 

Let X ~ be a space of the weights of faces whose coordinates are given by 
X := (xji)jgz/mz, i&z/nz- Let p be an embedding map from C(A’) to C(A1) given by 

(HO]). 

We define three types of actions Ri and Ri on df, corresponding to the three 
directions that the honeycomb bipartite graph can be arranged into rows. For j = 
1,... ,m, define x{i) := (xj,*,a:j+i,j)jgz/nz- Let Rj be the action on X induced by 
Rx(j) flll.2p with L = n. More precisely, Rj(x) = x' is given by 

' + cji) l=j, 

^ Xj-ij{l + u;-\)-i l=j-l, 

x,+i,(l + a;ri)-i l = j + l, 

^xii otherwise. 














50 


REI INOUE, THOMAS LAM, AND PAVLO PYLYAVSKYY 





Figure 15. Honeycomb dimer model on a cylinder. 


with oji := Ui{xj^i,... ,Xj^n) flU-iP - In a similar way, for i = 1,. .. ,iV, define x{i) : = 
{xm-j,N-h Xm-j,N-i+i)jez/mn'z and let Ri be the action on df given by 

with L = mn'. For i = 1, ..., M, dehne x{i) := i+ 2 -i, 

and let Ri be the action on X given by Rx{i) flll.2p with L = ^. 
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Figure 16. (n, m,/i;)-dimer 


Recall that we also have three types of actions on M given by the elements sj {j = 
0,1,..., m — 1) of the extended symmetric group W, the elements Sj (i = 0,1,..., iV — 1) 
of the extended symmetric group W, and the snake path action Tg {s = 1,, M). 


Proposition 11.3. 

(i) The actions Sj and Si are compatible with the actions Rj and Rj respectively: for 
xE X, we have 



P o RM) = s* o p{x 

), j E Z/mZ, 

(11.4) 


P o Mx) = s* o p{x] 

1, i E 1.1X1,. 

(11.5) 

(ii) For xeX, 

the action Pi satisfies 




p o = p(x), 

i = 1 ,..., M. 

(11.6) 

iii) For xE X, 

the snake path action T* satisfies 



T* op{x) = p{x), 

i = 1 ,..., M. 



Proof, (i) To show fill .dp . it is enough to prove the j = 1 case. We write Ei for the 
energy E{P~^QiP\ P~''Q 2 P^). The operator si acts on Ai as si{q) = q'\ 



E, 



Ei-i 


and the other qji do not change. By definition, Ri{x) = x' is obtained as 

+ j = 1, 

, ^ Xm,i+k{l + j = m, 

|x2,(l+a;ri)-i j = 2, 

^Xji otherwise. 


(11.7) 
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where cuj := On the other hand, by direct compntation we obtain 

<il,i+lEi , , _ g2,i+l-E’i+l 


p{^i) = 


n n TT^ ’ 

s=l - ns=l <lps 


p{l + OJi) — 




Thus we get 


(l2,i 


( -In , (i2,i+lEi+l 

P + iXi)) = - -= s 

P [Xm,i+kK^ + ^i-l) ) = - ^ - = ^1 

p[x2,i{l+u:i ) ) = -;;— = Si 


1 I - = Si op(a:i,i), 


Q.m,i+k+l \ 

/ 


Si ® piXm,i+]^ 1 


q3,iE.. 


i+l 


q2,i+i 

q3,i 


= Si op{x 2 ,i), 


and flll.dp follows. 

We can prove flll.Sp in a similar manner, by replacing Qi with Qi, P with P, and so 
on. 

(ii) Again, it is enongh to prove the case of i = 1. For simplicity, we write L for and 
set Xj := Xj^ 2 -j, Qj ■= Qjp-j for j G Z/LZ. From flll.dp . it follows that 

L 


p{Xj) = 


Qj 

qj+i'' 


P 


)=1 

w'=l 


( 11 , 8 ) 


Dehne Uj = ojj{xi ,..., xl) for j = 1,..., L. Using the dehnition of ujj and fill.81) . we 
obtain 

' 1 + \ Qi ( Ui 


P 


Ui-I 


qj+1' 


p 


1 + Wi 


= 1 , 


(11.9) 


by the following calculation: 


1+Uj _ 1 ns=i Xg + Xi ns=i Xi-s 


IXj-l 


Xi-1 ns=l Xi-l-s 

1 _ 1 I s^n-1 qi_t 
P ^ ^t=0 Qi+i Qi 


E n Qi-t 

i=l Oi 


UJn 


Qi 

^L-l T-rt 


Qi+l 


El=0 n.=l Xi-s 


1 . 


1 + Uj 1 - X, + Xi Ef=o^ nl=i Xi-s 
Thus we have po i?,j(i)(x(l)) = p(f(l)), and flll.bp follows. 

(hi) The snake path action Tg changes qj^i and in the same way if i + j = F + j' 
mod M. This condition is satished by gyj+i and qj+i,i, then the change is cancelled in 
p{xj^i). Thus we see that Tg o p{xj^i) = p{xj^i). □ 


References 

[BK] A. Berenstein and D. Kazhdan, Geometric and unipotent crystals, Geom. Funct. Anal., Special 
Volume, Part I, 188-236 (2000). 

[Eti] P. Etingof, Geometric crystals and set-theoretical solutions to the quantum Yang-Baxter equation. 
Comm. Algebra, 31, no. 4, 1961-1973 (2003). 

[Fay] J. D. Fay, Theta functions on Riemann surfaces, Lecture notes in Math., 352 (Springer-Verlag, 
Berlin-New York, 1973). 

[FM] V. V. Fock and A. Marshakov, Loop groups. Clusters, Dimers and Integrable systems, 
arXiv:1401.1606. 






























GENERALIZED DISCRETE TODA LATTICES 


53 


[GSTV] M. Gekhtman, M. Shapiro, S. Tabachnikov, A. Vainshtein, Integrable cluster dynamics of 
directed networks and pentagram maps, arXiv:1406.1883. 

[GK] A. B. Goncharov and R. Kenyon, Dimers and cluster integrable systems, Ann. Sci. Ec. Norm. 
Super. (4) 46, no. 5, 747-813 (2013). 

[HTI] R. Hirota, S. Tsujimoto and T. Imai, Difference scheme of soliton equations, in Future directions 
of nonlinear dynamics in physical and biological systems (Lyngby, 1992), NATO Adv. Sci. Inst. 
Ser. B Phys., 312, 7-15 (Plenum Press, New York, 1993). 

[IKT] R. Inoue, A. Kuniba and T. Takagi, Integrable structure of box-ball systems: crystal, Bethe 
ansatz, ultradiscretization and tropical geometry, J. Phys. A 45, no. 7, 073001, 64 pp (2012). 

[IT] R. Inoue and T. Takenawa, Tropical spectral curves and integrable cellular automata, Int. Math. 
Res. Not. IMRN 27 Art ID. rnn019, 27 pp (2008); 

-, A tropical analogue of Fay’s trisecant identity and an ultradiscrete periodic Toda lattice, 

Gomm. Math. Phys., 289, 995-1021 (2009). 

[ILP] R. Inoue, T. Lam, and P. Pylyavskyy, in preparation. 

[108] S. Iwao, Solution of the generalized periodic discrete Toda equation, J. Phys. A 41, no. 11, 115201, 
15 pp (2008). 

[110] S. Iwao, Solution of the generalized periodic discrete Toda equation. 11. Theta function solution, 
J. Phys. A 43, no. 15, 155208, 9 pp (2010). 

[109] S. Iwao, Linearisation of the (M, iL)-reduced non-autonomous discrete periodic KP equation, 
arXiv:0912.3333. 

[Kir] A.N. Kirillov, Introduction to tropical combinatorics, in “Physics and Gombinatorics 2000” (Eds. 
A. N. Kirillov and N. Liskova), Proceedings of the Nagoya 2000 International Workshop, 82-150, 
World Scientific, 2001. 

[Kho] A. Khovanskii, Three Lectures on Newton Polyhedra (Newton polyhedrons and Singularities). 
RIMS Kokyuroku, 1233, 1-17 (2001). 

[KKMMNN] S.-J. Kang, M. Kashiwara, K.G. Misra, T. Miwa, T. Nakashima, and A. Nakayashiki, 
Affine crystals and vertex models, in Infinite analysis Part A (Kyoto 1991), 449-484, Adv. Ser. 
Math. Phys., 16 , World Sci. Publishing, River Edge, NJ, 1992. 

[KNl] M. Kashiwara, T. Nakashima, and M. Okado, Affine geometric crystals and limit of perfect 
crystals. (English summary) Trans. Amer. Math. Soc. 360, no. 7, 3645-3686 (2008). 

[KN2] M. Kashiwara, T. Nakashima, and M. Okado, Tropical R maps and affine geometric crystals. 
Represent. Theory 14, 446-509 (2010). 

[KNY] K. Kajiwara, M. Noumi, and Y. Yamada, Discrete Dynamical Systems with 1E(A^1_]^ x 
Symmetry, Lett. Math. Phys. 60, no. 3, 211-219 (2002). 

[KT] T. Kimijima and T. Tokihiro, Initial-value problem of the discrete periodic Toda equation and its 
ultradiscretization. Inverse Problems, textbfl8, no. 6, 1705-1732 (2002). 

[Las] A. Lascoux, Double crystal graphs. Studies in memory of Issai Schur (Ghevaleret/Rehovot, 2000), 
95-114, Progr. Math., 210, Birkhauser Boston, Boston, MA, 2003. 

[LPll] T. Lam and P. Pylyavskyy, Affine geometric crystals in unipotent loop groups, Representation 
Theory 15, 719-728 (2011). 

[LP12a] T. Lam and P. Pylyavskyy, Total positivity in loop groups, I: Whirls and curls. Adv. Math. 
230, no. 3, 1222-1271 (2012). 

[LP12b] T. Lam and P. Pylyavskyy, Inverse problem in electrical cylindrical networks, SIAM J. of 
Applied Math., 72, 767-788 (2012). 

[LP13] T. Lam and P. Pylyavskyy, Crystals and total positivity on orientable surfaces, Selecta Math. 
(N.S.) 19 , no. 1, 173-235 (2013). 

[MIT] J. Mada, M. Idzumi and T. Tokihiro, Path description of conserved quantities of generalized 
periodic box-ball systems, J. Math. Phys. 46, no. 2, 022701, 19 pp (2005). 

[Rho] B. Rhoades, Cyclic sieving, promotion, and representation theory, J. Combin. Theory Ser. A 117, 
no. 10, no. 1, 38-76 (2010). 

[Sch] M.-P. Schiitzenberger, Promotion des morphismes d’ensembles ordonnes, Discrete Math. no.2, 
73-94 (1972). 

[Spe] D. Speyer, Perfect matchings and the octahedron recurrence, J. Algebraic Combin. 25, no. 3, 
309-348 (2007). 

[TS] D. Takahashi and J. Satsuma, A soliton cellular automaton, J. Phys. Soc. Japan 59 , no. 10, 
3514-3519 (1990). 

[vMM] P. van Moerbeke and D. Mumford, The spectrum of difference operators and algebraic curves, 
Acta Math. 143 (1-2) 94-154 (1979). 



54 


REI INOUE, THOMAS LAM, AND PAVLO PYLYAVSKYY 


Rei Inoue, Department of Mathematics and Informatics, Faculty of Science, Chiba 
University, Chiba 263-8522, Japan. 

E-mail address: reiiyOmatli. s . chiba-u. ac . jp 

Thomas Lam, Department of Mathematics, University of Michigan, Ann Arbor, MI 
48109, USA. 

E-mail address: tfylam@uinich.edu 

Pavlo Pylyavskyy, School of Mathematics, University of Minnesota, Minneapolis, 
MN 55414, USA. 

E-mail address: ppylyavs@umn.edu 


